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ABSTRACT 

The  effect  upon  the  correlation  functions  of  a 
fluid  of  adding  a  pair  potential  of  arbitrary  range  to  a 
repulsive  hard  core  is  investigated  in  the  case  of  the 
Percus-Yevick  approximation  of  classical  equilibrium 
statistical  mechanics.   The  nonlinear  P.Y.  equation  is 
perturbed,  to  first-order,  and  the  resulting  linear  integral 
equation  for  the  first-order  correlation  functions  is  then 
solved  exactly  for  a  general  Interparticle  potential.   The 
method  of  solution  is  presented  for  one-dimensional  systems, 
but  three-dimensional  non-polar  and  polar  liquids  can  be 
encompassed.  Essential  to  the  technique  is  the  calculation 
of  a  certain  zeroth-order  quantity,  which  is  proportional 
to  the  Faltung  of  the  hard  core  pair  correlation  function  with 
itself,  and  which  is  composed  of  a  sum  of  nth-neighbor 
contributions.   Obtaining  this  quantity  allows  conversion 
of  the  first-order  perturbed  P.Y.  integral  equation  to  an 
ordinary  differential  equation  of  fourth  order  for  the  direct 
correlation  function  inside  the  core.   Reduction  of  the 
inhomogeneous  term  and  determination  of  the  constants  of 
Integration  by  application  of  the  continuity  conditions  at 
the  hard  core  boundary  results  in  explicit  expressions  for 
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the  first-order  direct  correlation  function  inside  the 
core,  and  the  first-order  pair  correlation  function 
outside  the  core.   The  validity  of  the  first-order  perturba- 
tion theory  is  tested  by  evaluating  the  expression  obtained 
for  the  first-order  direct  correlation  function  in  the 
particular  case  of  a  tall  potential  which  is  a  replica  of 
the  core  over  a  finite  range.    A  linear  dependence  results. 
In  agreement  with  exact  statistical  mechanics.   A  simple 
application  is  made  to  the  glue  model  of  Baxter.   If  the 
phenomenon  of  linearity  is  extrapolated  as  a  model  for 
arbitrary  tail  potential,  and  the  P.Y.  equation  is  further 
simplified  to  the  Debye-Huckel  model  with  hard-core  insertion, 
a  number  of  explicit  solutions  may  be  obtained.   This  is 
carried  out  for  Yukawa  tall,  double  Yukawa,  and  power  law 
tails. 
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I.   Introduction 

The  Percus-Yevick  approximate  Integral  equation  for 
the  pair  correlation  function  In  classical  statistical  mechanics 
has  been  solved  exactly  for  hard  cores  [13].   The  direct 
correlation  function  Is  given  In  Its  non-vanlshlng  domain  by 
a  low  order  polynomial,  of  degree  one  or  three,  for  the  case 
of  one  or  three  -  dimensional  hard  cores.   Exact  solutions 
of  the  P.Y.  equation  for  potentials  of  arbitrary  range  have 
never  been  obtained,  however.   The  object  of  this  work  Is  to 
examine  the  effect  upon  the  correlation  functions  of  adding 
a  long-range  Interaction  to  the  repulsive  core.   Although 
treated  exhaustively  in  one-dimension,  our  technique  is  by 
no  means  limited  to  this.  Indeed  being  applicable  to  both 
three-dimensional  non-polar  and  polar  liquids. 

We  commence,  in  Sections  II  and  III,  by  reviewing  the 
formalism  of  distribution  functions  in  the  equilibrium  theory 
of  classical  fluids  [1,6].   The  pair  correlation  function  is 
defined,  and  its  role  in  determining  the  thermodynamics  of  an 
N-particle  system  is  discussed.   By  utilizing  a  functional 
Taylor  expansion  combined  with  the  relation  between  the 
density  of  a  nonuniform  system  and  the  pair  correlation 
function  of  a  uniform  system,  approximate  Integral  equations 
for  the  pair  correlation  function  can  be  derived  [3].   By  this 
technique,  the  Percus-Yevick  equation  is  obtained.   In 
Section  IV,  the  pair  potential  is  written  as  the  sum  of  two 
parts,  a  repulsive  hard  core  plus  a  tall  potential  of  arbitrary 
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range.   Making  use  of  the  fact  that  the  P.Y.  equation  has 
been  solved  exactly  for  hard  cores,  we  perturb  to  obtain  the 
first-order  Percus-Yevick  equation.   This  is  followed,  in 
Section  V,  by  a  detailed  description  of  the  method  of 
solution  of  this  first-order  perturbed  integral  equation. 
Expressions  for  the  first-order  direct  correlation  function 
and  pair  correlation  function  are  given.   In  the  following 
section,  the  perturbation  theory  is  tested  by  evaluating 
the  expression  obtained  from  it  previously  for  the  correlation 
function  in  the  case  of  a  hard  core  tail  potential  of  finite 
range.   Comparison  of  the  first-order  correlation  function 
with  its  counterpart  from  exact  statistical  mechanics  is  made 
and  the  results  are  discussed.   The  last  section  contains  a 
number  of  explicit  solutions  for  potentials  consisting  of  a 
hard-core  plus  tail.   First  order  phase  transitions  are  found 
to  occur  for  each  model  studied. 


-2- 


II.   Distribution  Functions  in  Classical  Equilibrium 
Statistical  Mechanics. 

We  begin  by  reviewing  the  formalism  of  distribution 
functions  in  the  theory  of  classical  fluids.   For  the  present, 
we  consider  a  single-component  simple  fluid  whose  total 
Internal  potential  energy  is  given  by  the  sum  of  pair 
interactions.   Later,  the  method  of  distribution  functions 
will  be  extended  to  a  more  general  type  of  potential  energy. 
The  assumption  of  two-body  interparticle  forces,  viz.,  that 
the  force  between  any  pair  of  particles  is  Independent  of 
the  presence  of  any  additional  particles,  has  as  a  consequence 
that  the  equation  of  state  depends  only  on  the  two-body 
distribution  function  and  not  upon  higher  order  distribution 
functions . 

A.     Coordinate  Distribution  Functions  in  a  Grand  Canonical 
Ensemble . 

Inasmuch  as  it  will  be  necessary  to  vary  the  density 
arbitrarily,  use  of  the  grand  canonical  ensemble  is  indicated. 
Thus,  if  the  fraction  of  the  systems  of  the  ensemble  having  the 
composition  N  in  the  region  d  r  d  £_  of  phase  space  appropriate 
to  N  is  denoted  by :  pjj'^  (r  ,c.)  d  r  d  p  ,    then  the  fraction 
of  systems  in  the  ensemble  with  composition  N,  viz.,  the 
probability  that  the  open  system  does  contain  exactly  N 
particles,  is  merely: 
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Pj"   (tjE.)  d  £  d  p.   The  integration  is  over  all  N-canonical 
coordinates  r,  .  .  .  r „  and  momenta  2_-,  .  .  .  2.m>    ^^'^   ^^    independent 
of  time.   Then 

(2.1)  I       [  ?jf\r,^)    d\  d%  =  1  . 
N=0  ■' 

(G ) 
The  density  function  ?},        can  be  written  as 

(2.2)  Pii^^r,p)  =   ^  e    ^  "  "  X^  , 

where  3  =  (kT)~  ,  and  X  is  the  fugaclty,  related  to  the 
chemical  potential,  y,  through   X=  e^  . 
H„  is  the  N-body  system  Hamiltonian, 

N  p?     ,    N  N 

1=1         lpj=l  1=1 

in  which  the  terms  are,  respectively,  the  kinetic  energy, 
internal  potential  energy,  and  external  potential  energy. 
The  normalizing  factor  Z^  is  the  grand  partition  function, 
given  by : 

(2.4)  Z      =      I  e  Xdrdp 

^   N=0  J 

(2.^a)        =  I      X^  Q 
N=0      ^^ 

where  Q^  is  the  canonical  ensemble  partition  function, 

(2.5)  Q^  =  -ofj^ —   e  d  r  d  p  . 
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The  momentum  integrations  in  (2.5)  are  trivial.   It  may  be 
verified  that 


(2.6)  Q,,   = 


^N 


N    "    N!    a3N    ' 


with   A   = 


h 


/27TmkT 


The  quantity  Z  ,  the  configuration  integral,  is  then  an 
N-fold  coordinate  space  integral. 


(2.7) 


^N  = 


e-''-^-^'--^-^r,...dr,  , 


in  which  V^  is  the  potential  energy  portion  of  (2.3) 
Introducing  the  activity,  z,  defined  by 


(2.8) 


z  = 


we  may  then  write  (2.4a)  as 


(2.9) 


N 


In  this  form  each  N-particle  system  of  the  ensemble  is  weighted 
by  z^/N! 

We  next  define  distribution  functions  expressing  the 
probabilities  of  observing  different  configurations  of  sets 
of  s  particles  out  of  the  total  number  N  (N  >_  s).   Let 
n^,^"^  (r-i  .  .  .r  )  dr-,...dr   be  the  probability  that,  if  the 

i^        _L       S        -1-         o 

configuration  of  the  system  of  N  particles  is  observed,  a 
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particle  (not  necessarily  particle  1)  will  be  found  in 

dr.,  at  r-,  ,  a  second  in  dr„  at  r,,,  ...,  and  another  in  dr   at  r 

—1         —1  '  —d  —d  b  u. 

Then,  in  the  canonical  ensemble. 


(2.10)   njj^^r^...r^)  =  j^ 


e  dr^  +  i-.-dllN 


^N 


with  Zj^  defined  by  (2.7) 


For  a  grand  ensemble,  we  denote  by  n  (r-,  .  .  .r  )dr-j^ 
the  probability  of  observing  particles  in  dr-,,...,dr   at 


.dr 


-1 ' • • • ' -s 


n  (rT...r  )  is  averaged  over  N  and  is  thus 
s  — 1   — s  ° 


independent  of  N.   Specifically, 


(2.11)   n  (r  .  ..r  )  =  I     n^^^r  ...r  ) 

si      S      j^^g    1^J       X      S 


Pj$°)(r,p)d\d%  . 


The  function  n  (r-,...r  )  is  the  s-body  distribution  function 


in  the  grand  canonical  ensemble;  it  depends  also  on  T  and  y 

tion  to  r, , . . • ,r  . 
— 1 '    '  — s 

From  (2.2),  (2.5),  (2.6)  and  (2.8),  it  follows  that: 


in  addition  to  rT,...,r 
— 1 '    '  — s 


z'^Z 


(2.12) 


N 


N!Z 


G 


Pji   (£>P)  d  r  d  p 


Hence,  using  (2.10)  and  (2.12),  (2.11)  may  be  written  as: 


(2.13)   n„(r,...r„)  =  Z"""  I     .^ 

N=s 


S  — 1    — S 


r  -3V^(^i---I1n^ 


dr3+i..drj^. 
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The  normalization  equation  for  n  (r-^^...r^)  Is 


(2.14) 


n 


(Li»- • -Lg)  ^^i* • -^^s 


V  N! 


pf)(r,p)  a\a\ 


\  (N-s)!  -^ 


B. 


The  Pair  Correlation  Function, 


The  dlmenslonless  s-body  correlation  function 

g  (r,  .....r^)  Is  defined  by 
°s  — 1     '— s 


(2.15)   ng(r^,. . .,r^)  =n^(r^)-"  n^ (r^ )gg (r^ , . . . ^r^ )  . 

The  above  equations  are  valid  for  both  crystals  and  fluids; 
we  now  specialize  to  the  latter.  For  a  uniform  system,  the 
one-body  density,  n-,  (r-,  )  is  a  constant. 


(2.16) 


'^l^'-l^  =  n  = 


_llm 


_N_ 
V 


=  const 


From  (2.15)  and  (2.16),  the  two-body  distribution, 
n„(r-,,rp)  becomes:   n2(r-|^,r2)  =  n  g2(ll]_jE.2^'   Because  of 
the  Isotropy  (spatial  translation  -  invarlance)  of  the  fluid, 
the  function  Ep^^i  >^p)    depends  only  on  r  =  |rp-r-|^|,  and  so: 
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(2.17)  n^Cr^jT^)  =  n^g(|r^-r2l)  . 

The  quantity  g (  I  H-i -£.p  I  )  is  known  as  the  pair  correlation 
function,  or  radial  distribution  function. 
Inasmuch  as  g(r)  determines  the  thermodynamics,  it 
is  the  principal  quantity  desired  in  the  N-body  problem 
of  classical  equilibrium  statistical  mechanics.   From 
(2.1^),  g(r)  is  asymptotic  to  unity,  in  the  thermodynamic 
limit,  except  during  a  phase  transition  when  N  has  a  double 
peak. 

The  total  energy,  E,  of  a  fluid  of  N-particles 
interacting  palrwise  with  a  potential  (J)(r)  in  the  absence 
of  external  fields  can  be  shown  to  be  [1] 

00 

(2.18)  E  =  I  NkT  +  I  n  I  (j)(r)  g(r)  ^iur^   dr  . 

0 

The  pressure,  P,  is  given  in  terms  of  the  grand  partition 
function: 


(2.19)  BVP  =  £n  Z_ 


In  expressing  P  directly  in  terms  of  g(r) ,  however,  there 
exists  both  the  virial  pressure  relation, 

00 

(2.20)     PV_  =  1  _  ^  I  r(!,.(r)  g(r)  ^i^r^   dr  , 

0 
and  the  compressibility  formula  [2], 
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(2.21) 


n  , 
6  ^T 


=  1  +  n 


[g(r)-l]  h-nr      dr  =  [1-  n 


2    -^ 
c(r)  hiw   dr]   , 


In  which  k   is  the  isothermal  compressibility ,  given  by 


(2.22) 


k  =  h^) 


The  quantity  c(r)  appearing  in  (2.21)  is  an  additional 
correlation  function  to  be  discussed  in  Section  111. 

The  pressures  calculated  from  (2.20)  and  (2.21)  will 
be  identical  if  the  exact  g(r)  is  used.   Otherwise,  there 
will  be  a  discrepancy  between  them  which  is  a  measure  of 
the  accuracy  of  the  theory  of  the  approximate  g(r)  used. 

Especially  appropriate  in  a  first-order  theory  of 
the  correlation  function  (e.g.,  see  equations  (4.4)  and  (4.11) 
below)  is  the  approximation  to  the  free  energy  F,  written  out 
through  second  order. 


(2.23)  F    =    Pq    +    I  nN 


CO 

r 


(t)-j^(r)    gQ(r)    47Tr      dr 


+    If  nN 


cl)^(r)    g-|_(r)    4iTr      dr    , 


where  F„  is  the  hard  core  free  energy.   The  energy  E  and 
pressure  P  then  follow  from  the  relations 


(2.24) 

and 

(2.25) 


^  =  h^^ 


P  =  _  ^  =   i 
3V   ~  B 


H^e 
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III.    Integral  Equations  for  the  Pair  Correlation  Function. 

We  next  proceed  to  examine  the  formal  procedures  by 
means  of  which  integral  equations  for  the  radial  distribution 
function  are  obtained  [31.   These  equations  provide  an 
approximate  theory  of  g(r)  for  the  N-body  problem. 


A. 


Functional  Taylor  Expansion. 


(3.1) 


First,  consider  a  regular  functional  [4]  of  degree  N 


N 


F[y(t)]  =  k^  +   ^ 


n=l 


b   b 
r  t 


a   a 


k^(t^...t^)y(t^)...y(t^) 


dt,  .  .  . dt 


where  the  kernel  k  (t,...t  )  is  symmetric  with  respect  to  its 
n  variables  in  the  interval  from  a  to  b.   This  can  be  extended 
to  an  infinite  series  by  letting  N  ->  <»,  such  that  F[y(t)] 

will  be  convergent  and  will  correspond  to  a  definite  value 

12        1/2 
y  (t)  dt]     is  less 

than  some  radius  of  convergence. 


(3.2)   F[y(t)]  =  k^  +  I 

^        n=l 


b 
r 


k^(t^. ..t^)y(t^)...y(t^)dt^..dt^ 


a 


The  variational  derivative  of  (3.2)  is 


5F 


(3.3)  j^-   k,(t,)  -^   I   n 

J   1'  n=2    •' 


k  (t.  , 
n   1 


t^)y(t„)..y(t„) 


dtp ... dt  , 
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Continuing  to  differentiate  p-1  times. 


«Pf 


(3-«)  6y(tp):..6y(t,)  'P'  ^(h---'p> 


b   b 


n=D+l         i         i 


a   a 


.dt 


Hence , 


(3.5) 


6^F 


6y(t  )...6y(t^) 


=  p!k  (t^...t  )  , 


y=o 


the  relation  between  the  kernel  k  (t-j^.-.t  ),  and  the  functional 
F[y(t)]  which  generates  it. 

The  functional  Taylor  expansion  can  then  be  written 


00 

(3.6)   F[y(t)]  =  F[0]  +  I      iy 


p  =  l 


fiPF 


a 


6y(t  )..  .6y(t^) 


y=o 


•  y(t^). . .y(tp)  dt^. . .dtp 


If  y  is  a  function  of  a  vector   r,  (3.6)  then  obtains  for 
F[y(r)].   The  formula  (3.6)  allows  us  to  obtain  expressions 
for  the  distribution  functions  and  thermodynamic  variables  in 
nonuniform  systems,  about  their  uniform  system  values  [5]. 
The  non-uniformity  can  arise  either  from  an  applied  external 
potential,  or  from  a  change  in  the  density;  F[y(r)]  is  then 
expanded  In  terms  of  whichever  function  is  chosen  to  be  varied. 
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B. 


The  Generating  Functional. 


We  now  recall  [6]  that  Z_  is  the  generating  functional 
for  the  complete  set  of  s-body  joint  distributions. 


^s^^l' 


,r  )j  (s  =  1,2,3,...)  given  in  (2.13).   First,  it 


-s 


will  be  convenient  to  employ  the  notation  (from  (2.9)  and  (2.3)) 


00     N 

(3.7)     Z  (0)  =  I     ^ 

^       N=0  ^^• 


-  I   I   Mr.,r  ) 


dr,  , 


,dr 


■N 


and 


(3.8) 


N 


00    N  r 


'G 


-2    .1      *(^i'^i) 


i?^J=l 


N 


e  dr.,  .  .  .  dr jT 


with  the  dimensionless  y(£.-)  "  ~  3U(r.  ),  to  denote  the  grand 
partition  function  in  the  absence  and  presence,  respectively, 
of  any  applied  external  potential,  U.   Similarly,  from  (2.13), 


3 


N 


N 


(3.9)   n^(r^...r^)  =  z;:ToT  J^  iT^   . 


-  f   I      Mr.,r  ) 


dUs+i-'-driN 


and 


3 


N 


(3.10)  n^(r^...r3|Y)  =  zTTVT  J 


00      N 

z 


gWT  ^i,  N-s! 


'   -n..l/'-'^^> 


N 

I  yCr.) 


From  (3.8) , 
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6 


N 


6Z, 


(3.11) 


SyCr-j^) 


N 

N^l  N-1' 
N 


I   YCr.) 
1=1    ^ 


dr 2 . . . dv^    , 


so  that 


6e 


N 


N 


,   -I  .J  ^^(r_.,r    )    M  yCr^) 


l?^j  =  l 


1=1 


dr2...dr^ 


Differentiating  a  total  of  s-tlmes,  one  finds 


6^Z  -Y(r.)      -Y(r  )   -    N 

(3.13)   7 ^ — 7 r  =  e   ~^  ...  e    ^    I        ^ 


fie 


.  .  .  6e 

N 

I 

i7^j  =  l 


N=s 


N-s! 


I  ,J    J^-^l'Zj^       Jj^^i^ 


N 

I 
1=1 


drg+i-.-dr^ 


-Y(rT)      -Y(r  ) 
=  e     ^  ...  e     ^ 


Zq(y)  n^(r^...rjY) 


Prom   formula    (3.13),    we   have,    for   s    =    1: 


(3.14)      n^(rjY)    =    z-t^e 


-^  Y(ri^)         5Z 


G 


6e 


f^    V     =    ^^-4 V    ^n    Z^(y) 

yCr^        6y(E.i) 


6n    (r    Iy)  S^'z    (y) 

^^^^'         fiYCrg)      =    Z^T^      SY(r2)fiY(r^)    "   "i (^^ I Y)ni(r2 I y) •    But 
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<5^Z^(y) 


from    (3.11),    z^TYT  6y(r^)6y(r^) 
Thus, 


=    n2(r^,r2lY)+   n^  (r^  |  y)  <5  (r^-r2  ) 


6n, (r, Iy) 
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In  general   [7] , 


(3.16)    — ^ 


6n    (r    ,. . . ,r    Iy) 

6Y(r)  q+1   —1'         '-q'— '' 


+   n^(r^,...,r^|Y)    I    \      ^i^^-^)    "   n-L(rlY)] 


i=l 


From  (3.7)  and  (3.8),  it  is  useful  to  note  the  following 


3 


N 


00        N  00        N 

(3.17)    Z    (0)    =      I      It  2^(0)    =      I      ^ 

^  N=0       •       '^  N=0    ^  • 


-   f        I      *(r.,r    ) 


3 


dr^...dr^ 

N 


00        N_s  00  ]\i_s 

y     ^ z       (0)  =     y     ^ 

,/      N-s!    ^N-s^^^  ./  N-s! 

N=s  N=s 


^    ^  i7^J=S  +  l  ^       <^ 


dUs+l-.-dr^ 


and 


00  N  <»        N 

(3.18)   z   (Y)   =     I     ^r^N^Y)  -     I     W  , 

^  N=0  •       ^^  N=0    ^^-     J 

N 


N 


r      -f   .J      *^^1'^.1^ 


i?^j  =  l 


J 


I   Y(r.) 

i=l  ^  ^ 

e  dr    .  .  .  dr j^ 
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N 


°o    N-s 


N=s 


N-s!   N-s 


«>    N-s  f 


N=s 


N-s! 


N 


I   Y(r,) 


For  s  =  1,  (3.9)  and  (3.10)  yield 


N 


N 


(3.19)  n,(rj  = 


1-1'  -  Z^TOT  j^i^  N-1! 


r      -t  .J   *^^i'^.l) 


i?^j  =  l 


dr. 


'^^N 


N 


(3.20)  n^(r^lY) 


1     y   ^ 

Z^TTy^il  N-1!  J 


.  ^-ij^^^^^^^p 


N 

I  Y(r.) 


1  =  1      ^      ^ 
e         dr „ . . . dr„ 


It  is  advantageous  to  relabel  as  follows. 


N 


,     0°    N-1 
(3.21)  n^(r2)  =  z-^T  ,J,  hzY 


,     -  I  .J  J^^^>L.) 


±9^j=2 


dr-. 


'^^N 


B 


N 


00 

(3.22)  n^CrjlY)  =  ^4-y  I 


N 
.I^Y(ri) 

•  e         dr  T  .  .  .  drpr  , 


which  can  then  be  interpreted  as  the  one-body  density  in  an 

(N-l)-particle  system  with  an  external  potential  energy, 

N 

I      y(r.).      Similarly, 
i=2     ^ 
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(3.23)  n^(rg) 


0°    N-s  +  1  r 


V^N=s   N-^  = 


drg+l-'-drj^ 


and 


N 


(S-ZI)  n,(rjy)  =  j^^f-T  J 


N-S  +  1 
1     V   z 


G'^'  N=s 


N-s! 


N 


.1  yCr.) 


1  =  S 


drg^l-'.drj^ 


are  the  one-body  densities  in  an  (N-s+l)-particle  system 

without  and  with,  respectively,  an  external  potential  energy, 
N 

I  Y(rii). 
i=s 

Furthermore,  the  non-uniform  two-body  distribution 

can  be  written 


B 


N 


(3.25)  n2(r2,r3|Y)  =  ^4^  J 


N-1  r 
z 


G^Y^  N=3  ^-3' 

N 


e^        dr^...drj^ 
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c. 


Bootstrap  Operation 


It  is  now  possible  to  show  [8]  that  there  exists  a 
relation  between  uniform  and  non-uniform  distribution  functions, 
provided  that  the  applied  external  potential,  Yj  is  of  a 
special  type,  denoted  by  Y'  ^^   particular,  yCe.,-)  is  the 

J 

external  potential  created  at  r.  by  fixing  s-particles  at 
the  positions  Z-\  >  •  '  '  >L     ^^'^   allowing  the  remaining  (N-s)- 
particles  to  interact  pairwise  with  the  s-centers  of  force 
each  sees ,  i.e.. 


(3.26)    Y(r^,) 


I      *(ll.>£i) 


J=l 


is  the  external  potential  energy  of  a  particle  at  r . .   For 

J 

instance,  from  (3.24) 


,      °°  N-s 

(3.27)  n-L(r^^jY(r3^l))  =  Z^TyT  ^J  +  i  N-(s  +  l)! 


N 


N 


-      2       .        J  ^^Hl'Hn-)     .1    Y^^i) 


i/j=s+l 


i  =  s  +  l 


dr 


■s  +  2 


•dr^  > 


where,  from  (3.l8) 


00    N-s 
(3.28)  Z(y)    =       I       ^ 


G 


N=s 


N-s! 


. 

3 

r        -     2 

•    •    • 

e 

N  N   ^ 

i^j=s+l    ^   J  ^i=s+l 


drg+i-.-dr^ 
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The  quantity  n,(r  +i Iy^Hq+i^)  ^^  ^^^  one-body  density  in  an 
(N-s )-particle  system,  with  an  external  potential  y  created 
by  fixing  particles  at  r-i  >•••)£  . 

An  important  special  case  occurs  when  only  a  single 
particle  is  fixed,  at  r -^ .   Then,  putting  s  =  1  in  (3 .  26)-(3  .  28) 
gives 

(3.29)  n-^(r2lY(r2)) 


00  N-1 

y  2 

i.  N-2 ' 

N=2  ^  ^• 


N 


N 


,      -I  .J   *(lli>£j)   .1  Y(r,) 


i?^j  =  2 


1  =  2 


dr- 


^^N 


«>   N-1  r 

N^i  ^-1' 


N 


N 


,   -  I  .  I   <|.(r.,r  )   J  YCr^) 


i?^j=2 


1  =  2 


dr„, 


'^^N 


with 


(3.30) 


Y(r^)  =  -  B(l)(r^,r^) 


The  function  n-,  (rp  1 7(^2) )  is  the  one-body  density  in  an 
(N-l)-particle  system,  with  an  external  potential  y   created 
by  fixing  one  particle  at  r, .   Because  of  the  interaction 
between  the  fixed  particle  and  the  remaining  particles,  the 
density  about  the  fixed  particle  is  nonuniform.   In  this  sense, 
the  distribution  of  matter  about  a' fixed  narticle  can  be  thought 
of  as  the  response  of  a  system  to  an  "external'  field,  y. 
Also,  from  (3.25) 
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e 


N 


(3.31)  n2*^^2»-3'^^ 


N-1 


h^  nk   ^-3'  i 


N 


I   YCr.) 


with  y(£.-)  given  by  (3.30). 

Now,  in  addition  to  the  ability  to  introduce  the 
potential  y  1^^  which  fixed  particles  play  a  special  role, 
another  consequence  of  the  procedure  of  setting  aside  fixed 
coordinate  locations  r-,,...,r   is  that  conditional  probability 

J.  o 

distributions  are  introduced  into  the  interpretation.   That 
is,  "fixing"  particles  is  a  formal  device  which  allows 
conditional  distribution  functions  to  be  introduced  [9]. 
For  instance,  the  one-body  density  n-,  (rp  |  y  (^2) )  is  merely 
the  probability  density,  P2(ll2llli)'  °^  finding  a  particle 
at  rp  when  it  is  known  that  there  is  a  particle  at  r -,  . 
That  is,  by  fixing  a  particle  at  r-,  ,  the  one-body  distribution 
n-,  (rp  I yCjlp) )  becomes  a  conditional  two-body  distribution 


(3.32) 


n^(r2|Y(i:2))  =  P2(^2'-l'* 


These  conditional  distributions  provide  the  connection 
between  uniform  distributions  and  non-uniform  distributions, 
viz.,  with  Y  turned  on.   Thus,  from  (3.9) 
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(3.33) 


1^2^-1 '-2^ 


n 


I^ 


N 


N-1 

y  2 

L       M_2 1 

N=2  '^  ^- 


-I  -J  .'^^-I'^-P 


l/j  =  l 


dr-, .  .  .dr 


•N 


N 


N-1 

y  ^ 

N=l   ^^  ^- 


,   -|.J_*(r,,r) 


i^j  =  l 


dr^. 


,dr 


•N 


having  multiplied  both  numerator  and  denominator  by  z 


.  N                .    N  N 

-  I  .1  /(-i>rj)  -  I     ^*(ri,rj)   -  3.1  *(i:i>-i) 

D  ^.  1?^J=1                        l7^j=2          ''  1=2 

But  e  "^          =  e      "^          -e 


Fix  one  particle  at  r , .   Then,  from  (3.30) 


I 


N  =  2 


N-2! 


N 

I 
l7^j=2 


»   _N-1  ,    ,  -I  ^J_J(L,>LP      ,.L  Y(r,) 


N 

I 
1  =  2 


dr. 


•^^N 


N 


N 


N^l  N-1' 


-  f  .  ?   Mr.,r  )   J   Y(r.) 


l?^j=2 


1  =  2 


drp, 


'^^N 


Direct  comparison  with  (3.29)  reveals  that 
(3.35)  \^IJ        =n^(r2lY(r2))  . 


A  corresponding  relation  may  be  obtained  for  s-centers  of 
force  at  r-,,...,r  .   Thus,  again  from  (3.9) 
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(3-36)  n^ir_^,...,._^, 


N-s 


N 


J+1    N-(s  +  l)!    . 


,      -I       I      Krii,r    ) 

— s  +  2 


N 


00  N-s 

V  z 


N=s 


N-s! 


.         -f    .J       *(^1>^,-) 


1/J=1 


dr    , -, 
— s+1 


■'^^N 


having  multiplied  both  numerator  and  denominator  by  z 
But  e    1^J=1      '   =  e    ^^J-^1      ' 

N      S  R     ^ 

•  e   ^  ^^^  J  ^        .  e    ^'^J  .  Fix  s-particles 


at  £.1  J  •  .  .  jT  .   Then,  from  (3.26) 


(3.37) 


■^s  +  l^-l'  •  •  •  '-s+1^ 


^s^^l"--'^s^ 


N 


N 


N-s 


N=Ll  N-(s  +  l)!  . 


,   -  f    I        *(ri,r  )    I   Y(r  ) 


^^s+2- 


•^^N 


00     N-s 

,/    N-s! 
N=s 


N 

I 
i7^j=s+l 


,        -   I  .^.I^^^C^i.^j)  .  L  jCr.) 


N 
i=s+l 


^^s+l-'-'^^N 


I  .  I  _*(!:,>£,) 


i?^1  =  l   '^   "^ 
The  factor  e      "^  ,  which  represents  the  pairwise 


interaction  of  the  s-centers  of  force  among  themselves,  has 
cancelled  out. 
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Comparing  the  right  hand  side  of  (3.37)  and  (3.27) 
indicates  that 


(3-38)       °;\  ',...,/)'   =  "l(Es+ll^(!ls+l>' 

o    -L  o 


The  one-body  distribution  with  s-particles  fixed  at  r  , .  .  .  ,r 
is  a  conditional  (s+l)-body  distribution,  p  . .  (r  ,-,|r,  ,...,r  ) 

(3.39)    n^(r3^J9(r3+i))   =  P^^^(r^^^\r^,-^-,r^)    • 

In  the  first  order  approximate  theory  of  the  pair  correlation 
function,  it  is  necessary  to  fix  only  one  particle;  thus  (3.35) 
is  the  equation  which  will  be  utilized. 

The  significance  of  (3.38)  and  (3.35)  is  that  they 
enable  higher-order  distribution  functions  to  be  evaluated 
in  terms  of  lower-order  ones  when  the  potential  y   is  turned  on; 
that  is,  they  offer  a  "bootstrap"  operation  [10]. 

At  this  point  we  remark  that  the  external  potential  y 
employed  in  deducing  formula  (3.13)  could  have  been  y ;    the 
same  result  would  follow  with  suitable  labelling  to  designate 
the  locations  of  the  fixed  centers  of  force.   In  particular, 
the  relation  (3.15),  in  the  presence  of  a  single  center  of 
force  at  r, ,  could  be  written 


6n-|^(r2lY(il2)) 


(3.^0)   =  np(rp,r^|Y)  -  n-,  (r,  |  Y)n-,  (r^  |  y) 

6y{v^)  ^  -^  -3       1-2    1  -3 


+  n-^(r2lY)6(r2-r2) 
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in  which  (refer  to  (3.1^)) 


(3.41) 


n-j^{r^\y{r^))    = 


6Zq(y) 


Zq(y)   5y(I12) 
1     ^^G^^^ 


Z^(y)   6y  (L3) 


with 


(3.42)   Y(rL2^  =  -  ^^(Li>Z2^     >        ^^£3)  =  "  3  (|)  (  T;^  ,  Tg  ) 


D 


Integral  Equations 


1.    The  Direct  Correlation  Function. 

The  formulas  (3.6),  (3.35),  and  (3.4o)  enable  various 
integral  equations  for  the  pair  correlation  function  to  be 
obtained.   First,  suppose  the  function  to  be  varied  is  the 
applied  external  potential,  and  that  the  one-body  density  is 
expanded  as  y  is  turned  on.   Then,  from  (3.6) 


(3.43)     n^(r|Y)  =  n  + 


6n-j^(r  1y) 


6Y(iL) 


Y(y_)  d;^  + 


Y=0 


Turn  on  a  particle  at  x,  for  example.   Then  for  fluids, 
Y(r)  =  -B())(r-x),  y{^)    =   -3<j)(y-x)   a  function  of  relative 
coordinate   alone. 

Equation  (3.35)  then  implies 


(3.44)     n^(r|Y) 


n2(r-x) 
n 


n^{y-x) 


n^CylY)  =     n 


by  virtue  of  (2.l6).   Next,  (3.4o)  can  be  used  to  write 
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<5n,  (r  Iy) 

=  n„(r,y|Y)  -  n.,  (r  |  Y)n-,  (y  |  y)  +  n,  (r  I  y)  6  (r-y ) .  Thus 

<SY(iL)  i  -    i  -       i 

(3.^3)  becomes 


n„(r-x) 
(3.^5)   -^ =  n  - 


2 

[n^(r-y_)    -  n  +  nfi  (r-y  )  ]  3())  (y-x)  dy  + 


This  as  a  Debye-Huckel  equation  [11],  valid  for  weak  potentials, 
(|)(r).   Upon  Fourier  transforming,  (3.^+5)  may  be  solved  in 
the  form 


(3.46)        (g-l)j^  =  = 


1  +  nB<j)j^ 

Alternatively,  if  the  density  is  varied  instead  of  the  external 

potential,  a  simpler  indication  of  the  degree  of  inhomogenelty 

will  be  given.   It  is  thus  necessary  to  evaluate  the  "inverse" 

6Y(r-,) 

derivative,  — = 7 1 — v  :  the  effect  of  a  density  change  on 

'   <5n-|_(r2lY) 

any  quantity  in  the  system  may  then  be  found  by  means  of 

the  chain  rule.   Inspection  of  (3.15)  suggests  that  a-      , 1 — r- 

should  also  be  split  into  singular  and  non-local  parts: 

<SY(r;^)  1 

(3-^7)   5^^(^^|^)   =   n^(rglY)   ^^^1-^2^  "  ^^^l'^2\'<^    ' 


This  then  expresses  6Y/6n-,  in  terms  of  the  dimensionless 
quantity  C,  which  is  called  the  direct  correlation  function, 
From  the  chain  rule. 
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(3.48) 


-SCZi-Ilg) 


Substituting  (3.15)  and  (3.^7)  Into  (3.^8)  and  dropping  the 
explicit  specification   |y)   yields 


(3.^9)   g(lli>Il2^  "^  "  '^^Zi>'L2^    =   n(r2)c(r-j_,r2)[g(r2»^3^~^^^^3 


Equation  (3.'^9)  is  an  exact  Interrelation  between  g  and  c, 
When  iterated,  (3.^9)  becomes 


(3.50)  gir-^,r^)-l   =    c(r-j^,r2)  + 


c(r^,r2)n(r2)c(r2,r2)  dr^ 


c(lli3E.3)i^(l^3)c(ro)Ili|)n(r^)c(r|^,r2)  dr^  dr|^  +  •• 


The  right  hand  side  of  (3-50)  is  the  deviation  of  g  from 
its  non-interaction  value  of  unity. 

2.    The  Percus-Yevlck  Equation 

For  the  case  of  rapidly  varying  or  discontinuous  potentials 
—  the  extreme  case  would  be  a  rigid  sphere  interaction  with 
an  infinite  discontinuity  —  then  the  functional  Taylor  expansion 
of  n-,(r|Y)e~^  —  with  respect  to  n-,(r|Y))  when  truncated  at 
the  linear  term,  will  be  quite  accurate  due  to  the  effect  of 
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the  smallness  of  the  direct  correlations  in  the  higher  terms 
Thus,  from  (3.6) 


(3.51)   n^(r|Y)e"'»'^^^  =  n  + 


6n^(r|Y)e"^^-^ 


<Sn3^(y|Y) 


[n-|_(y  I  Y)-n]ciy, 


n-^^(y)=n 
(y  =  0) 
Using  (3.^'+)  and  (3.^7),  and  choosing  the  origin  as  x,  yields 


(3.52)   g(r)e 


3(}>( 


^)-l  =  n  I 


c(r-r')(g(r' )  -1)  dr 


However,  for  uniform  systems  (y  =  0),  and  (3.^9)  becomes 


(3.53)   g(r)  -1  -  c(r)  =  n 


c(r-r')(g(r')  -1)  dr'  , 


which  is  a  general  relation  between  g  and  c  for  a  uniform 
system. 

Thus,  from  (3.52)  and  (3.53),  the  approximation  (3.51) 
can  be  expressed  in  the  form 


(3.5^^) 


c(r)  =  (1  -  e^'^^^h    g(r) 


This  is  the  Percus-Yevick  approximation.   Using  (3.53) 
and  (3.5^)  to  eliminate  c(r)  results  in  a  single  integral 
equation  for  g(r), 

(3.55)   g(r)e^*^^^  -  1  =  n  j  (g(r-r')  -1)(1  -  e^*^^'^)  g(rM  dr' 
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Equation  (3.55)  is  known  as  the  Percus-Yevlck  (P.Y.) 
equation  [12].   Defining  the  correlation  function 
T(r)  =  g(r)e  ''^  —  ,  and  using  the  notation  e(r)  =  exp  (.-^<\>ir))  , 
f(r)  =  e(r)  -1,  (3.55)  may  be  written 


(3.56)     T(r)  -1  =  n 


[T(r-r' )e(r-r ' )  -1]  T(r')  f(r')  dr'  , 


a  nonlinear  Faltung  equation. 

The  quantities  e(r)j  f(r),  g(r)  and  c(r)   have  finite 
discontinuities  at  any  separation  at  which  (i)(r)  is  discontinuous 
However,  T(r)  is  always  continuous.   The  relation  between 
Cj  gj  and  T  in  the  P.Y.  approximation  is,  from  (3.5^)  and  the 
definition  of  T(r), 

(3.57)         c(r)   =   g(r)  -  T(r)  . 

From  (3.5^),  it  is  noted  that  the  range  of  c(r)  is  exactly 
that  of  the  Intermolecular  potential  (})(r),  in  the  P.Y. 
approximation . 
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IV.    The  First-Order  Perturbed  P.Y.  Equation 

The  nonlinear  integral  equation  (3.56)  has  been 
solved  exactly  for  the  case  of  hard  rods  (one  dimension)  and 
hard  spheres  (three  dimensions)  [13].   That  is,  the  potential 
4)(r)  has  the  form 

(^.1)    4>  =  0   (|r|  >  a)  ,    (})  =  <»   (|rl  <  a)  , 

where  a  is  the  range  of  the  potential.   Denoting  this  type 
of  potential  by  the  symbol  <|)Q(r),  we  have  e„(r)=  exp(-3<))„  (r ) )  , 
fQ(r)  =  ^^(r)  -  1,  and  the  correlations  CQ(r),  Sn^-^  ^"^^  '^0^-^' 
related  in  the  P.Y.  approximation  by 

{H.2)        CQ(r)  =  TQ(r)fQ(r)  ,    gQ(r)  =  TQ(r)eQ(r)  . 

It  is  now  convenient  to  refer  to  (3.56),  with  the  potential 
(^.l)j  as  the  hard  core  P.Y.  equation,  or  zeroth-order  P.Y. 
equation.   Written  out  in  terms  of  the  zero  subscripts, 

(^.3)    TQ(r)  -  1   =   n(TQeQ  -  1)  *  r^f^  . 

The  pair  correlation  function  g„  can  be  expressed  as  a  sum 
of  nth  -  neighbor  contributions,   and  the  direct  correlation  Cq 
(the  negative  of  which  continues  g^  into  the   hard  core  interior) 
is  given  by  a  low  order  polynomial,  of  degree  one  or  three  for 
the  case  of  one  or  three  dimensional  hard  cores.   The  difference, 
gQ-  Cq,  gives  Tq.   The  explicit  form  of  these  functions  will  be 
written  down  in  the  subsequent  development,  when  required. 


-28- 


Exact  solutions  of  C3.56)  for  general  potentials  of 
arbitrary  range  have  never  been  found.   It  is  however,  possible 
to  make  progress  toward  this  goal  by  exploiting  the  fact  that 
the  solutions  of  the  zeroth-order  equation  (4.3)  are  known. 
The  appropriate  course  would  be  to  perturb  the  nonlinear 
equation  (3.56),  and  then  attempt  to  solve  exactly  the  resulting 
first-order  perturbed  equation. 

We  write  the  interaction  potential  cf)(r)  as  a  sum  of 
two  parts,  the  repulsive  hard  core  potential  (})„(r),  plus  a 
tail  potential  ())-,(r),  which  can  be  of  arbitrary  range: 

(4.4)  (l)(r)    =    (l)Q(r)    +    (()-^(r)    , 
where 

(4.5)  c})^   =    0      (|r|    <   a)    ,  <t>j_  ^   0      {\r\    >    a)    . 

Since  adding  an  Infinite  range  attractive  type  of  potential, 
(p-.  ,  would  certainly  be  expected  to  modify  the  effects  of  the 
short-range  core  potential,  4)^,,  it  is  of  dubious  value  to 
regard  <()-,  as  being  "small  enough"  to  be  a  perturbation  on  (j)q. 

The  quantity  occurring  directly  in  (3.56)  is  not  the 
potential  <f)(r),  but  its  exponential,  e(r).   This  suggests  the 
possibility  of  perturbing  (3-56)  in  powers  of  the  exponential 
of  ())-,,  rather  than  (j),  itself.   Indeed,  since  Q^i^)    is  a  unit 
step  function,  and  exp  (-3<J)-,(r))  approaches  unity  rapidly 
(particularly  for  (t)-,(r)'s  which  ^  0  rapidly  as  |r|  ->■  °°)  ,  the 
quantity  e-,  (r)  ,  where 
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(4.6)  e-^(r)  =  e(r)  -  e^ir)    =   eQ(r)[exp  (-&<p-^ir))    -1] 

can  properly  be  regarded  as  a  polntwise  perturbation  on  e^(r) 
for  either  attractive  or  repulsive  <(),  ,  if  -  e  <3(}'-|  <   —   • 
Then,  to  first-order  in  this  perturbation  scheme, 

(4.7)  T(r)  =  TQ(r)  +  T^(r) 

in  which  t-,  (r)  is  the  solution  of  the  first-order  perturbed 
P.Y.  equation: 

(4.8)  T-|_(r)  =  ndgeQ-  1)  «  (TQe^+  x^fg)  +  nd^e-j^t  T^e^)  *  x^f^. 

Equations  representing  successively  higher  orders  of  the 
perturbation  can  also  be  obtained;  the  ultimate  justification 
for  a  first-order  perturbation  will  be  a  test  of  the  closeness 
of  agreement  between  the  results  it  yields  and  the  exact 
results . 

As  opposed  to  the  nonlinear  (3.55)  on  (3.56),  (4.8)  is 
now  a  linear  equation  for  T-,(r),  which  may  be  cast  in  the 
form  of  a  Fredholm  Integral  equation 

(4.9)  T-j^(r)  =  n  j  k(r,r')  T^(r')  dr'  +  P(r) 

in  which  the  kernel  k(r,r')  is  given  by 

k(r,r')  =  TQ(r-r')[eQ(r')fQ(r-r')+  e^ (r-r ' )f Q(r ' ) ]-f q (r ' ) , 

(4.10)  and 

F(r)   =   n(ToeQ  -  1  +  r^fQ)  *  x^e^  . 
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Clearly  the  above  kernel  is  not  a  function  of  (r-r')  alone, 
meaning  that  (4.9)  Is  not  a  convolution  equation.   Consequently, 
immediate  solution  of  (4.8)  by  transforming  and  inverting  is 
unfortunately  not  possible,  and  the  task  of  solving  becomes 
non-trivial. 

The  P.Y.  equation,  perturbed  to  first-order,  can 
alternatively  be  expressed  in  terms  of  c-,(r)  and  g-|(r),  rather 
than  T-,(r),  as  is  (4.8).   To  do  this,  first  write,  as  in  (4.7) 

(4.11)  g(r)  =  gQ(r)  +  g-^(r)  ,     c(r)  =  CQ(r)  +  c^(r)  . 
Then,  from  (3.53), 

(4.12)  g^(r)  -  c^(r)   =   n  Cq  *  g^  +  n  c^  *  (gQ  -1) 
and  from  (3.54), 

(4.13)  c^(r)  =  gQ(r)e-^(r)  ,     (|r|  >  a)  . 

Equations  (4.12)  and  (4.13)  are  an  expression  of  the  first- 
order  perturbed  P.Y.  equation. 
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V.      Solution  of  First-Order  Perturbed  P.Y.  Equation: 
One  Dimension. 

A.    Conversion  to  a  Differential  Equation. 

A  considerable  effort  was  devoted  to  obtaining  a 
technique  for  the  exact  solution  of  the  first-order 
perturbed  P.Y.  equation  for  general  potentials  of  arbitrary- 
range   [1^4].   Each  method  devised,  although  interesting  in 
Itself,  was  invariably  discovered  to  be  deficient  in  some 
significant  way.   We  shall  not  discuss  these  shortcomings 
here,  but  Instead  proceed  to  explain  the  successful  technique 
finally  developed.   This  latter  technique  is  suitable  for 
both  non-polar  and  polar  liquids. 

1.   Sth  -  Neighbor  Structure  of  y{x)    and  r(x). 

We  begin  by  observing  that  equation  (^.12)  can  be 
written  in  the  form 

(5.1)  g3_(r)  =  r  *  C^ 
where 

(5.2)  r(r)  =  (6(r)  +  n(gQ-l))  *  (5(r)  +  n(gQ-l))  , 

Itself  a  Faltung,  is  yet  to  be  calculated. 

Equation  (5.1)  Is  verified  by  Fourier  transforming  (4,12), 

(5.3)  g^(k)  -  c-j^(k)  =  n  CQ(k)g-j^(k)  +  n  c^(k)  (gQ-1)  (k)  . 
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Rearranging, 


1   +   nCin-DCk)      . 


>1^-^ 


1   -  nc^Ck) 


But    from    (3.53), 


-1 


(5.5) 


1   +   n(g„-l)(k)    =    (1   -   nc„(k))       . 


Thus,  (5.^)  may  be  written 


(5.6) 

having  denoted 

(5.7) 


^-^(k)  =  r(k)  c^(k)  , 


-2 


r(k)  =  (l-ncQ(k))    =  (l+n(gQ-l)(k)) 


and  (5.1)  follows  by  transforming  (5.6)  back  to  coordinate  space 
It  is  of  interest  to  calculate  first  y(^)j  defined  as 


(5.8) 


Y(r)  =  6(r)  +  n(g  -l)(r) 


In  one  dimension,  the  Fourier  transform  pair  for  y(x)  and  yik)  is 


Y(x)  =  2^   [  e-^^^  Y(k)  dk 


(5.9) 


Y(k)  = 
From  (5.8)  and  (5.5), 


ikX    /   ^   -, 

e     y(x)  dx 


(5.10) 


y(x)   =   ^ 

2tt   J 


-ikx 


1  -ncQ(k) 


dk  . 


For  one  dimensional  hard  cores,  the  explicit  expression  for  the 
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symmetric  function  Cq(x)  is  found  [13]  from  solution  of 


(^.3)  to  be 


c„(x)  =  -  Q^(l-nlx|)  ,     |x|  <  a 


(5.11) 

=0  ,     |x|  >  a 

where   Qq  =  (l-na)~  . 

For  such  a  one  dimensional  hard  core  fluid,  the  exact 
c„(x)  can  be  shown  [15]  to  have  precisely  the  range  a,  so 
that  (5.11)  is  not  an  approximation  but  is  identical  with 
the  exact  result. 


Using  (5.11),  the  quantity  l-nc„(k)  is  readily  calculated 


to  be 


(5.12)     1  -  ne^(k)    =   A^dc)  +  A^(k)  e"^^^  +  A_(k)e^^^  , 

in  which 

2(nQ.)^ 
A,(k)  =  1  +      ° 


nQ       (nQ„)^ 
(5.13)  A^(k)  =  ^  i 5— 

nQ      (nQ^)^ 
A  (k)  =  -  -f-^  i  -     ° 


k  2     ' 


SO  that,  from  (5.10) 

00 

-ikx 


(5.14)    y(x)  =  ^ 


00 

e 


2Tr 


_l   AQ(k)+  A^(k)e-^^^+  A_(k)e^^^ 


dk  . 


Since  (Aq+  A^e  ^'^^+  A_e^^^)  is  non-vanishing  (see  (5.21))  and 
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a  continuous  periodic  function  of  period  2tt  in  the  variable 
(ka) J  a  Fourier  series  expansion  can  be  performed 


-1     °°   ^ 
(5.15)   (A„Ck)  +  A. Ck)e-^^^  +  A  (k)e^^^)    =    I   G  (k)e 


-iska 


where  the  Fourier  coefficients  G  (k)  are  determined  from 


27r 


(5.16)   Gg(k)  =  27 


iska 


Q  AQ(k)  +  A^(k)e-^^^  +  A_(k)e^^^ 


d(ka)  . 


ika   ^ 

With  the  substitution  z  =  e    .  G  (k)  can  be  written  as  a 

s 

contour  integral  around  the  unit  circle. 


(5.17)   G^(k)  =  2¥r-  ; 


=  1  '^^o'^  ^-^"^  *+^  '^^ 


dz 
z 


Suppose    s    >    0.      Then 


(5.18)      G    (k)    =   ^ 

S  2lTl 


^^    (A_z      +    AqZ    +    A^) 


dz 


For  s  <  0,  it  is  necessary  only  to  complex  conjugate  G  (k) , 
that  is. 


(5.19) 


G   (k)  =  G  (k)  . 


The  poles  of  the  integrand  in  (5.18)  occur  when  z 
z  =  Zp,  where 


(5.20) 


z,  and 


^1  = 

2A_ 
-    2A      - 

Aq-4a_A^ 
2A_ 

^2    = 

/Aq-4a_A^ 
2A 
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Inserting  the  values  of  Aq,  A_,  and  A^  from  (5.13)  into  the 
square  root  reveals  that 


(5.21) 


/Aq  -  ^A_A^  =  1  . 


Further  use  of  (5.13)  implies  that 


(5.22)    z^  = 


k 


1  +  i 


"2  =  1  -  ^  nQ 


nQ 


0 


0 


One  pole  (that  at  z-,  )  is  inside  the  unit  circle,  whereas 

the  second  pole  is  outside.   Therefore,   G  (k)  =  2tti  Res  (z,) 
■^  s  i 


There  follows 


(5.23)    G^(k)  =  z^ 


And  from  (5.19) , 


s 


1  +  i 


nQ, 


(s  >  0) 


k 


(5. at)   G^(k)  =  (1  -  1  j^) 


From  (5.1^)  and  (5.15), 


(s  <  0) 


(5.25)   y(x)  =   I     ^ 


S  =  _oo 


2-n 


f  ^-ik(x+sa)  -  (j^)  ^^    ^ 


or 


(5.26) 


y(x)  =    I     Gs(x+sa)  , 


where  G  (x)  =  0  for  s/x  >  0  (refer  to  discussion  for  r"(x) 
following  (5.37)). 

Thus,  the  well-known  conclusion  [13]  that  y{x),    and  hence 
&,q{x)  ,    is  expressible  as  a  sum  of  s   -  neighbor  contributions. 
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Returning  to  equation  (5.1)  and  the  quantity 

r(x)  =  Y  *  y(x)>  it  will  now  be  shown  that  r(x)  also  has 

an  S^"  -  neighbor  structure.   Since  rCk)  =  CyC^))  ,  we 
have,  corresponding  to  (5.10), 


(5.27)  r(x)    = 


2Tr 


-Ikx 


(1   -  ncQ(k))' 


dk 


In   view   of    (5.12) 


(5.28)  r(x)    = 


2tt 


00 

r 


i 

_oo 


-ikx 


(AQ(k)    +   A^(k)e-^^^   +   A_(k)e^^^)^ 


dk 


Analogous  to  (5.15),  there  follows 


(5.29)     (AQ(k)+  A^(k)e"^^^+  A_(k)e^^^)  ^=    I      rg(k)  e"""^^^ 


_2 

I 


in  which  the  Fourier  coefficients  r  (k)  are  given  by 


(5.30)    r^(k)  =  27 


2tt 


0 


iska 

ry— ry— — ^  d(ka)   , 

(AQ(k)+  A^(k)e"^^^+  A_(k)e'-^^)^ 


which  may  be  written  as  a  contour  integral  (choosing  s  >_  0 , 
for  example) 


(5.31) 


r  (k)  =  ^ 

S         27T1 


S  +  1 


^^  (A_z^  +  AqZ  +  A^)^ 


dz 


The  (double)  poles  of  the  integrand  of  (5.31)  occur  again  at 
z^  and  Z2,  given  by  (5.22).   Residue  evaluation  of  (5.31) 
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results  in 

2  2 
n  Q  k  "^ 

(5.32)  r  (k)  =  Cs  +  1  +  2  — ^)C1  +  i  -3-)    ,   (s  >  0) 

For  s  £  0,  complex  congugate  r  (k) .   Thus, 

2  2 

n  Q  ,   s 

(5.33)  r  (k)  =  (-  s  +  1  +  2  — ^)(1  -  i  -§-)    ,   (s  <  0) 


From  (5.28)  and  (5.29)  , 

oo 

(5.34)    r(x)  =  ^   1   f  e-i^(^+sa)  p  ^^^  ^j^ 


2Tr   ^ 

S  =  _oo        ' 

—  00 


S 


or 


(5.35)    r(x)   =    ^  r  (x+sa)  ,   where  r  (x)  =  0  for  |  >  0, 


S  =  -oo 


a  sum  of  s^'^  -  neighbor  contributions,  which  will  now  be 
examined  in  greater  detail. 


2.    Evaluation  of  r"(x). 


From  (5.32)  and  (5-33),  it  is  seen  that  in  addition 
to  the  poles  of  order  s  at  k  =  nQ„i  (for  s  >_  0)  and  k  =  -nQ^i 
(for  s  <_  0),  there  exists  a  double  pole  at  k  =  0 .   This  latter 
pole  is  undesirable  from  several  points  of  view,  and  will  be 
avoided  by  considering  instead  r  (x+sa).   Thus,  from  (5.35), 

o 
00 

(5.36)    r"(x)  =   I      r"(x+sa)  ,       (_oo  <  x  <  °°) , 

3  =  _oo 
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where 


(5.37)    r^(x+.a)  =  -  ji  I  e-i«>=+==')(k2  ?^Ck))  dk  . 

—  00 

such  that  (5.32)  is  to  be  used  for  s  >  0,  and  (5.33)  for  s  <  0 
To  evaluate  (5.37),  the  semicircular  contour  must  be  closed 
in  the  upper  half  plane  if  (x+sa)  <  0 ,  and  closed  in  the  lower 
half  plane  if  (x+sa)  >  0.  Since  the  algebraic  sign  of  s  fixes 
the  location  of  the  pole,  there  will  be  an  extended  range  of  x 
for  which  a  given  r^(x+sa)  vanishes,  becoming  and  remaining 
non-zero  outside  this  range.   In  concise  form. 


(5.38)    r^(x+sa)  =  0  ,    s(x+sa)  >  0  , 


-°°  <  X  <  oo 


r"(x)  is  a  symmetric  function  which  is  composed  of  an  infinite 
number  of  distinct  sth  _  neighbor  functions,  r"(x+sa).   For 
sa  <  |x|  <  (s+l)a,  r"(x)  is  a  finite  sum  of  (s+l)-terms; 
this  piecewise  character  is  a  computational  asset. 

The  first  few  r^(x+sa)  involve  only  delta-functions 
and  their  derivatives,  and  can  be  written  down  immediately 
from  (5.37), 

T'^ix)    =    -  2n^Q^    6(x)  +  5"(x) 
(5.39)     r;^(x)  =  -  2n2Q2  ^(x)  -  2nQ,  6'(x) 


'0  "^"^  ~  ^'^'^0 
Iq    6(x)  +  2nQQ 


r_;;(x)  =  -  2n2Q2  5(x)  +  2nQ^  5'(x) 
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3.    Correlation  Function  Inside  the  Hard  Cores 

It  is  now  possible  to  convert  equation  (5.1)  into 

a  differential  equation.   Thus,  suppose  x  lies  within  the 

range  of  the  potential  (0  <  x  <  a).   Since  g-,  (x)  vanishes 

in  this  range,  (5.1)  can  be  written 

a 
(5.^0)      I  r(x-x')  bj_(x')  dx'  =  -  A(x)      (0  <  x  <  a), 
-a 

where,  from  (4.13) 


(5.41)     A(x)   = 


r(x-x')  gQ(x')  e^(x')  dx' 


b-|_(x)  =  c^(x)   (|x|  <  a)  ,    b^ix)    =    0    (|x|  >  a) 

Once  b-j^(x)  is  determined,  it  is  then  possible  to  obtain  g-,(x) 
(=  r*b-j^+A(x))  ,  X  >  a. 

Taking  two  derivatives  of  (5.40),  and  noting  that 
the  argument  of  r"(x-x')  in  the  b^-integral  has  the  range 
-a  <  x-x '  <  2a,  in  which  range  r"  is  composed  solely  of 
delta-functions  (see  (5.39)),  leads  to   the  deduction  that 


(5.42)   b![(x)  -  2nQQb|(a-x)  -  2n^QQ[b^(x)  +  b^(a-x)]  =  -A"(x), 


a  second-order,  first-difference  equation  for  b.(x),  in  which 
the  primes  denote  derivatives  with  respect  to  the  arguments. 
By  letting  x  ^  a-x,  another  second-order,  first-difference 
equation  is  obtained  for  b,(a-x) 
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(5.^3)   b"Ca-x)  -  2nQQb|(x)  -  2n^QQ[b-^(x)  +  b-^(a-x)]  =  -A"(a-x). 


(5.42)  and  (5.43)  are  simultaneous  differential  equations 
for  b-,(x)  and  b-,(a-x).   With  D  =  g^  j  these  are 


(D^-2n^QQ)b^(x)      +  (2nQQD-2n^QQ)b^(a-x)=  -  D^A(x) 
(5.44a,b) 

(-2nQQD-2n^QQ)b^(x)  +  (D^-2n^QQ)b-^ (a-x)    =  -  D^A(a-x) 

Eliminating  b-,  (a-x)  results  in  a  single  differential  equation 
of  fourth  order  for  b-,(x),  valid  in  the  interval  0  <  x  <  a: 

(5.45)     D^b^(x)  =  (2n^QQA"  -  a'-^)(x)  -  2(n^QQA"  +nQQA '"  )  (a-x)  . 


B.    Solution  of  Differential  Equation. 

1.    Reduction  of  Inhomogeneous  Term. 

Denoting  the  right-hand  side  of  (5.45)  by  F(x), 
there  follows 

(5.46)    b^(x)  =  b^(0)  +  b^(0)x  +  |  b!|(0)x^  +  |  b^"(0)x^ 

X     X 

+  I  ••  •  I  F(x)  d  X  . 


0    0 
F(x)  and  each  of  the  four  constants  of  integration 

t  It  Kt 

b  (0),  b-,(0),  b-,(0),  and  b^  (0)  depend  on  the  tail  potential, 

(()-,(x).   Since  (5.46)  is  valid  for  all  x  in  the  interval  0  <  x  <  a 

these  four  constants  must  be  determined  identically  for  a  fixed 
'i)l(x). 
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The  process  of  reducing   the  multiple  Integral  in 
(5.46)  generates  another  cubic  polynomial  in  x.   When 
this  is  combined  with  the  cubic  already  appearing  in  (5.46) 
(before  reduction  of  the  multiple  integral).  It  then  follows 
that 


(5.47)   b-L(x)  =  (2n^QQ-D^) 


X  X 


0  0 


A(x)  d  X 


a-x  X 


-  (2n^QQ-2nQQD) 


A(x)  d  X 


+  Bq  +  B,x  +  B2X^  +  B-,x^  , 


(0  <  X  <  a) , 


and 


(5.48)   b^(a-x)  =  (2n2Q^-D^) 


a-x   X 


0    0 


A(x)  d  X 


-  (2n^QQ+2nQQD) 


X   X 
f  r 


a  a 


A(x)  d  X 


+  Bq  +  B^(a-x)  +  B2(a-x)^+  B^(a-x)-^,  (0<x<a), 


where  the  constants  B  ,  B-,  ,  Bp   and  B^  have  resulted  from 
the  reduction  of  (5.46),  and  must  be  determined  identically 


The  formula   •••   F(x)  dx...dx   =  -? ^-n- 

k        K  , ;    (n-1)! 


{x-O^-^   F(0  d? 


n  times     n  times 
can  be  useful  in  this  regard. 
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for  all  X  in  the  interval  0  <  x  <  a.   However,  before 
applying  the  boundary  conditions  on  b-,  (x)  to  accomplish 
this,  it  is  necessary  to  first  determine  whether  the  four 
B.  (j =0,1, 2, 3)  are  linearly  independent.   That  is,  although 
(5.^5)  is  of  fourth  order,  it  resulted  from  two  second-order, 
first-difference  equations  either  of  which  implies  the 
other  when  the  symmetry  of  b-,  (x)  is  used;  the  fourth-order 
differential  equation  is  equivalent  to  a  second-order,  first- 
difference  equation  for  symmetric  b-,(x). 

2.    Linear  Dependence  of  Integration  Constants. 

Thus,  first  substituting  (5.^7)  and  (5.^8)  into  (5.44a), 
and  then  into  (5.44b),  and  simplifying,  yields  respectively: 


(5.49)  [-2n 


4^4 


(a-2x)A(x)  dx  -  2n 


3^3 


A(x)  dx 


+   ^2   ~   n^o^^l^  ^^^2  +3a^B  )  -  n^QQ(2BQ+aB^+a^B2+  a^B^)] 


+  x[3B2+  2nQQ(B2+3aB2)+  n^QQ(2aB2+  Sa^B^)] 


and 


-  nQ^x  [3B  +  nQQ(2B2+  3aB  ) ]  =  0  ,    (0  <  x  <  a) 


(5.50)  [-2n^QQ 


a  a 

(a-2x)A(x)  dx  -  2n-^Q;^    A(x)  dx 
0  0 


+  (B2+  3aB2)-  hQqB^-  n^QQ(2BQ+  aB^+  a^B^^-  a^B  )  ]  + 
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+  xE-SB^-  2nQQB2+  n^Q^(2aB2+  3a^B^)] 


-  nQ^x^[3B^+   nQQ(2B2+  SaB^)]  =  0  ,     (0  <  x  <  a) 


Equating  to  zero  the  coefficient  of  each  power  of  x  in  (5.^9) 
and  (5.50)  furnishes  six  equations  for  the  four  B.'s.   Only 
two  of  these  equations  are  independent,  not  being  implied  by 
the  others.   Thus,  there  results  a  system  of  defect  two, 
given  by: 

(5.51)  -2n^QQBQ-  nQQ(l+anQQ)B^+  (l-a^n^Q^ )B2+  a{3-a.^n^Ql)^y   A 

2nQQB2   +    3(l+anQQ)B2   =  0 


where 

a  a 


A  =  2n^Qjj  I  (a-2x)  A(x)  dx  +  2n^Q,^ 


'0 


A(x)  dx 


0  0 


Assigning  B,  and  B-,  arbitrarily,  B^  and  Bp  are  then  linear 

combinations  of  B,  and  B-,,  given  by 

3(l+anQQ) 
^2  = 2HQ^  ^3 

^0  " T~2 2nQ     ^1 5~2^2Hi~^^"^^'^'^0^ 


(5.52) 


■(l-a^n2QQ)-a(3-a^n^QQ)]B2, 
B-,  ,  B,,  arbitrary. 


1'   3* 

(5.^7)  thus  contains  two  arbitrary  constants,  to  be  determined 
from  the  continuity  conditions  at  the  hard  core  boundary. 
It  is  written  finally  as 
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X 


a-x 


(5.53)    b^(x)    =    2n^QQ 


-    2nQ 


(x-c)A(?)    dC    -    2n 


2^2 


(a-x-c)A(c)    d? 


0 
a-x 


A(C)    dC    -    A(x)    + 


a 


+    B      +    B-^x    +    BgX^    +    B^x^    , 


(0   <   X   <   a) 


(5.53)  is  the  exact  solution  of  (5.1)  for  0  <  x  <  a. 


3.    Continuity  Conditions  at  Hard-Core  Boundary. 

It  can  be  shown  directly  from  (4.8)  that  in  one 
dimension  x-,  (x)  and  its  first  derivative  are  continuous  at 
X  =  a.   In  three  dimensions,  t.,  and  its  first  three  derivatives 
are  continuous  at  the  hard  core  boundary.   As  a  consequence. 


(5.54a,b) 


g-L(a  )  =  gQe-^^(a  )  -  b-j^(a  )  , 

'   +  '   +      '   — 

^^{b.    )    =    (gQe-]_)  (a  )  -  b-^(a  )  . 


Eliminating  the  g^  (a  )  and  g-,  (a  )  by  combining  (5.54a,b)  with 
the  basic  equation  (5.1)  yields  two  equations,  implicitly. 


for  B   B_.   These  are 


(5.55a,b) 


r(a'^-x')b^(x')  dx'  +  A(a'^)  =  gQe^(a'^)  -  b^(a  ;B^,B2) 
r'(a'^-x')b-^(x')  dx'+  A'(a"^)  =  (g^e^)  '  (a"^)-b-[_(a";B^  jB^) . 


-a 


Because  of  (5.55),  it  is  now  necessary  to  obtain  r(x)  itself; 
thus  far,  only  r"(x)  has  been  evaluated  (refer  to  (5-37)). 
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Knowledge  of  rCx)  is  also  a  requirement  for  calculating  g-,Cx) 
once  b-,Cx)  has  been  obtained. 

As  shown  in  the  Appendix, 


(A. 5)   r(x)  =  5(x)  +  2ncQ(x)  -  Sn^Q^ | x |  +  K  ,   (-2a  <  x  <  2a), 

with 

(A. 9)   K  =  an^QQ(6  -  6an  +  ^  a^n^  -  a^n^)   . 

Inserting  r(a''"-x')  and  r'(a''"-x')  from  (A. 5)  into  (5.55a,b), 
and  explicitly  reducing  (T*  b,)(a  )  and  (f*  b,)(a  )  leads  to 


(5.56a)   2[(K-3an^QQ-nQQ) 


I  b^(x',B^ 
0 


,B,)  dx' 


2^2 


x'b^(x', 6^,62)  dx']  =  gQe^(a'^)-  A(a^)-  b-^(a  ,6^,62) 


and 


(5.56b)   -^n^Q^  I  b^(x',B-^,B2)  dx'  =  (gQe^)'(a"^) 


-A'(a  )  -  b-^(a",B^,B2) 


We  next  use  (5.53)  to  calculate  b-,(a~,B  ,B_)  and  b-,(a  ,B  B^), 
and  then  substitute  these  into  (5.56).   The  result  is 


(5.57a)   2[(K-3an^QQ-nQQ)  [  b^(x)  dx  -  n^Q^ 


xb-,(x)  dx] 


Q 


o  2„2 
=  -2n  Qq 


a 

r 


(a-2x)ACx)  dx  -  2nQQ  |  A(x)  dx  - 
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-  ;^(a,nQQ)  -  [a  +  k(a,nQQ)]B^  -  [a^H-  3a  k(a,nQQ)+  m(a,nQQ)  iB^ , 


and 


2^2 


a 
r 


(5.57b)  -^^Qq        b^Cx)  dx  =  -2nQQ  A(0)  -  B^-  [3a  +6ak(a,nQQ)]B2 


In  which 


k(a,nQQ)  =  - 


(l+anQ^) 
2nQ„ 


A 


(5.58)    £(a,nQQ)  =  -  ^  2^2 


2n"Q 


0 


m(a,nQ„)  = 5-^  [5^  (1+anQ  )  ( 1-a^n^Q^ ) 

-   a(3-a2n2Q^)] 

With  the  help  of  (5.53),  the  integrals  on  the  left-hand  sides 
of  (5.57a)  and  (5.57b)  can  be  written 


(5.59) 


a 
r 


b-^(x)  dx  =  a^(a,nQQ)  - 


A(x)  dx 


and 
(5.60) 


+  2nQ 


0 


xA(x)  dx 


0 


a  a 

+  n^Q^Ca^    A(x)  dx  -  2a 


X  A(x)  dx] 


0 


a 

e 


X  b-|^(x)  dx  =  a2(a,nQQ)  - 


0 


X  A(x)  dx 


+  nQ, 


(2ax-x  )A(x)  dx 


+  2n2Q2  [|_ 


A(x)  dx  -  |—  I  X  A(x)  dx  - 


a 
2 


x^A(x)  dx  +  ^   x^A(x)  dx]  ; 


3  J 


0 
_2|7_ 


where 

2 

(5.61)  cj-^(a,nQQ)  =  a^(a,nQQ)  +  [|-  +  ak(a,nQQ)  ]B-^ 

+  [^  +  a^k(a,nQQ}  +  amCa,nQQ)]B2 

3    2 

(5.62)  a^(a,nQQ)    =  |  a^Jl(a,nQQ)  +  [|-  +  |-  k(a,nQQ)]B^ 

+  [|-  +  I  a^k(a,nQQ)  +  j-   m(a,nQQ)]B2 


^.    Explicit  Equations  for  the  Two  Arbitrary  Constants 

Substituting  (5.59)  and  (5.60)  into  (5.57a)  and 
rearranging  enables  us  to  write 


(5.63)         h^^B^  +  h^2B3  =  ^1 

with  the  notation 

(5.64a)  h^-^  =   2a(K-3an^QQ  -  nQQ)(|+  k(a,nQQ)) 

-  2a^n^QQ(|  +  |  k(a,nQQ))  +  a  +  k(a,nQQ) 

3 

(5.6ilb)   h-,^2  =  2a(K-3an^QQ  -  nQQ)(^+  a\(a,nQQ)  +m(a,nQQ)) 

-  2a^n^QQ(|  m(a,nQQ)  +  |-  a^k(a,nQQ)  +  |-) 
+  (a^+  3a^k(a,nQQ)  +  m(a,nQQ)) 


(5.64c)   R^  =  ,|^^  -  £(a,nQQ)[2a(K-  3an^QQ-  nQ^)  -a^n^Q^  +1] 
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(5.64d)  ^^   =    [(2K-6an^QQ-2nQQ)(l-a^n^QQ)  -  2nQQ 


2an  Qg  +  j   a^n  Qq]  •    ACx)  dx 

0 


2^2 


2^2 


,2^2^„„.3n3 


+  [(2K-6an''QQ-2nQQ)(2an^QQ-2nQQ)+  2n  QQ+ilan-'QQ 


-  2a  n  QqJ 


X  A(x)  dx  -(2n3Q^+2an  Qq) 


X  A(x)  dx 


+  3  n  Qq  J 


X   A(x)  dx  . 


0 


Substituting  (5.59)  Into  (5.57b)  and  rearranging  leads  to 


(5.65) 


^^21^1  "*"  ^^22^3   ~   ^^2  ' 


where 


(5.66a)   h 


21 


1  -  2a^n^QQ  -  4an^QQk(a,nQQ) 


(5.66b)   h 


22 


3a(a+2k(a,nQQ))  -  4a^n^QQ  (|  +  k(a,nQQ)) 


2  2 
-  4an  QQm(a,nQQ) 


(5.66c)   R„  =   i>o  -   2nQ„A(0)  +  4a  n^Q^il(a,nQ„ ) 


(5.66d)   ijj. 


-  ^-'^0 


a 

A(x)  dx  +  8n^Q^ 


+  4n\J  [a^ 


a 

r 


A(x)  dx  -  2a 


X 

A(x) 

dx 

0 

a 

X 

A(x) 

dx] 
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Thus,  we  have  manipulated  (5.55a,b)  into  two  simultaneous 
equations  —  refer  to  (5.63)  and  (5.65)  —  for  B^  and  B^, 
explicitly . 

From  (5.63)  and  (5.65),  B   and  B^  are  found  to  be 

h„„R-,   —  h-,  qRq 
(5.67a)       B^  =  -^^-^ ^^-^- 

^^11*^22  ~  ^12^^21 
(5.67b)       B,  =  'll'2  -^21% 


3    h-^^h22  -  h^2h2i 


With  b,(x)  determined,  g-,(x)  follows  from 

(5.68)         g^(x)  =  r  *  b-j^(x)  +  A(x)  ,    (x  >  a) 
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VI.     Comparison  with  Exact  Results. 

To  examine  the  validity  of  the  first-order  perturbation 
theory,  we  can  avail  ourselves  of  the  fact  that  the  exact 
correlation  functions  are  known  for  a  one-dlmenslonal  system 
of  hard  cores  (see  remark  below  equation  (5.11)).   What  Is 
therefore  Indicated  Is  an  evaluation  of  b-,  (x)  ,  given  by 
equation  (5.53),  for  the  case  of  "extended  hard  cores".  I.e., 
for   (j)-]  (x)   Itself  a  hard  core.   Geometrically,  this  Is  a 
means  of  decreasing  the  excluded  length  by  increasing  the 
size  of  the  cores  while  maintaining  the  density  constant. 
Once   b-,  (x)   is  evaluated,  comparison  can  be  made  with  exact 
statistical  mechanics. 

A.     Extended  Hard-Core  Kernel 

Before  performing  the  calculation  of  b-,  (x)   from 
(5.53)  J  which  is  tedious,  we  present  the  following  simple 
considerations  for  the  purpose  of  obtaining  a  quick  estimate. 

Regarding  extended  hard  cores  as  a  perturbation  about 

hard  cores,  the  first-order  change  in  direct  correlation 

function,   6c(x),   is  given  by 

00 

(6.1)         6c(x)  =  I  A(x,x')  e^(x')  dx '  , 

in  which,  for  finite  core  extension  from  a  to  a+b , 

e-|^(x'  )=-l,    a<|x'|<  a+b  , 
and  the  kernel  A(x,x')   is  to  be  determined.   From.  (5.11), 
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(6.2)    6c(x)  =  (l-n|x| )[ 


e(a-|xj)     e(a+b-|x| ) 


(1-na)'^     (l-n(a+b)) 


]  , 


where  e(x)  is  the  unit  step  function. 


(6.3) 


e(x)  =  1  ,    X  >  0 
=  0  ,    X  <  0 


Thus  J  (6.1)  may  be  written 


(6.^1) 


(l_n|x|)[  ^i^zM) ei^^b-\x\)    3  ^ 


(1-na)^     (l-n(a+b))' 


a+b 


A(x,x')  dx'  - 


-a 

r 


A(x,x ' )  dx 


-(a+b) 


Since  A(x,x')  is  even  in  both  x  and  x',  it  follows  that 


(6.5)   - 


a+b 
f 


A(x,x')  dx'  =  I  (l-n|x|)[  ^l^zMl  _  ^(a+b-|x|)^  ^^ 


(1-na)^    (l-n(a+b)) 


from  which 


(6.6) 


a 

A(x,x')  dx'  -  ^  (l-n|x|)  ^^^"'^'  j,   =  constant. 


0 


(1-na)' 


Differentiating  (6.6)  with  respect  to  a  yields 

(6.7)    A(x,a)  -  n(l-n|x|)  ^1^=1^1   _  ^  ilzUhll   6(a-|x|)  =  0; 

(1-an)-^    ^  (l-na)"^ 

and  therefore  the  result  for  the  extended  hard-core  kernel 
A(x,x')  is 
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(6.8)   A(x,x')  =  n(l-n|x| )  - 


where   | x ' |  £  — • 


(Ix'|-Ix|)  .  1  5(|x'|-|x|) 


We  shall  defer  use  of  (6.8)  until  the  last  chapter, 
and  proceed  next  to  the  evaluation  of  the  first-order  P.Y. 
correlation  function  (5.53)  for  a  hard-core  tail  potential, 

B.     Calculation  of   A(x) 


We  choose   <})-,  (x)   to  be  a  hard  core  potential  over 
the  Interval   a  <  |x|  <  a+Aa  ,   with  (extended)  range 
Aa  <  a.   For  all  other  values  of  |x|,   (t)-,(x)  is  zero. 
Consider  the  interval  0  <  x  <  a,  over  which  b-,(x)  is  non- 
zero.  We  want  to  calculate  A(x),  defined  as 


(6.9) 


A(x)  =   I  r(x-x')  gQ(x')  e^(x')  dx'  . 


From  (6.9), 

a+Aa 


(6.10)   A(x)  =   I    [r(x-x')  +  r(x+x')]  gQ(x')  e^(x')  dx ' , 

a 

(0  <  X  <  a) ; 


e-,(x')  =  -  1   (a  <  x'  <  a+Aa),   and  is  zero  otherwise, 
accounting  for  the  truncated  limits  of  integration  in  (6.10) 
For  the  arguments  of  F  in  (6.10),  there  is  the  implication 
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-(a+Aa)  <  (x-x')  <  0,   and   a  <  (x+x')  <  2a+Aa  . 

Thus,  I x-x 'I  <  2a   for  all  x  in  the  interval  0  <  x  <  a, 
and  all  x'  in  the  range  of  ({>,  ;  whereas,  (x+x')  can  excee<i 
2a  when  x  >  (a-Aa).   However,  for  its  argument  exceeding 
2a,   r  is  given  by  an  expression  different  from  that  which 
obtains  when  its  argument  is  less  than  2a.   Consequently, 
A(x)  will  have  different  functional  forms  for  0  <  x  <  a-Aa, 
and  a-Aa  <  x  <  a.   Systemization  is  achieved  by  changing 
variables  in  (6.10),  yielding 


x-a 


(6.11)   A(x)  =     j    r(y)gQ(x-y)e-^(x-y)  dy 


x-a-Aa 
x+a+Aa 


J    r(y)  gQ(y-x)  e^(y-x)  dy  ,   (0<x<a-Aa), 


x+a 


and 

(6.12)   A(x)  = 


x-a 
r 


x-a-Aa 


r(y)  gQ(x-y)  e-j_(x-y)  dy 


2a 

f 


x+a+Aa 


r(y)gr,(y-x)e,  (y-x)dy  + 


r(y)gQ(y-x)e-^(y-x)dy,  (a-Aa<x<a) 


x+a 


2a 


The  y-limits  of  integration   in  (6.11)   and  (6.12)  determine 
which  expression  to  use  for  r(y),  i.e.,  either 

(A. 5)     r(y)  =  6(y)  +  2ncn(y)  +  ^  -  Sn^Q^lyl,  (-2a<y<2a)  , 


-0 


or 


-5^- 


3   2^2, 


(B.7)   r(y)  =  6(y)  +  2nCQ(y)  +  K  +  |  h'^Qq  (  |  y-2a| -2  |  y  |  ) 


-  ^nQ^e 


-nQQ(y-2a)     ^    ^  -nQQ(y-2a) 


+  n  QQ(y-2a)e 


,  (2a<y<3a) 


(as  shown  in  Appendices  A  and  B). 

Consider  first  the  Interval  0  <  x  <  a-Aa.   Since 

c  (y)  =  0  when  |y|  >  a,  the  limits  of  intej^ration  in 
x-a 

•^n^y  ^Sr|(^~y)  "^y   will  depend  on  whether  0  <  x  <  Aa, 
x-a-Aa 
or  Aa  <  X  <  a-Aa.   If  we  choose  Aa  <  a/2   (similar 

considerations  obtain  for  Aa  >  a/2),  there  follows 


(6.13) 


x-a 
r 


x-a-Aa 


x-a 
r 


CQ(y)gQ(x-y)dy,  0<x<Aa 


Co(y)go(^-y)  ^y  =  "* 


-a 


x-a 


c  (y)g^(x-y)dy,  Aa<x<a-Aa, 


x-a-Aa 


Consequently,  A(x)  will  be  a  continuous  function  of  x  for 
0  <  X  <  a,  but  will  have  different  functional  forms  for  each 
of  the  three  intervals  0  <  x  <  Aa,  Aa  <  x  <  a-Aa,  and 
a-Aa  <  x  <  a  (such  that  0  <  Aa  <  a/2).   We  denote  these 
three  expressions  by  A-,(x),  A2(x)  and  A-,(x),  respectively. 
From  (6.11),  we  find  (for  extended  hard  cores) 
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x-a 


x-a 


(6.14)   A-^(x)  =  -2n   J  CQ(y)gQ(x-y)  dy  -  K        ggCx-y)  dy 


-a 


-  3n  Qq 


x-a 
r 


J 
x-a-Aa 


x-a-Aa 

X 

x+a+Aa 


y  gg^x-y)  (iy-  K 


gQ(y-x)dy 


x+a 


2  2 


x+a+Aa 


x+a 


y  gg^y-x)  dy. 


0  <  X  <  Aa 


and 


(6.15)   A2(x)  =  -2n 


x-a 


x-a 


CQ(y)gQ(x-y)  dy  -  K 


gQ(x-y)dy 


x-a-Aa 


x-a-Aa 


x-a 


x+a+Aa 


o    2^2 
3n  Qq 


y  gg^x-y)  dy 


-  K 


gQ(y-x)dy 


X-a-Aa 


x+a 


x+a+Aa 


2  2 


y  gg^y-x)  dy  , 


Aa  <  X  <  a-Aa 


x+a 


Working  out  the  integrals,  substituting  their  values  into 
(6.14)  and  (6.15),  and  simplifying  leads  to  the  result 


(6.16) 
with 

(6.17) 
and 


A^(x)  =  C-^  +  C^x    , 


0  <  X  <  Aa 


(6anQQ  +  6Qq 


^^)    +  [-6n(a+Aa)Q2-6QQ+  ^]e~''  °  ^ 


C„  =   2nQ 
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(6.18)         ApCx)  =  B-,  +  B^x  ,  Aa  <  X  <  a-Aa 

with 

B^=    2[(l+2an)QQ+  2Qq-  ^]-  2[ (l+2n(a+Aa) )Qq+2Qq-  ^]e 

Bg  =  2nQQ(  1  -  e    "^   )  . 

From  equation  (6.12),  after  evaluating  the  Integrals  and 
rearranging  terms,  we  obtain 

nQ„(x-a)      -nQ„(x-a) 

(6.20)  A^(x)  =  A^  +  A2X  +  A^e         +  A^e 

-nQQ(x-a) 
+  Ar-xe         ,  a-Aa  <  x  <  a 

with 

(6.21)  A^    =    [2(l+2an)QQ    +    4Qq    -   ^] 

11  c:  PR  v:      -nQp.Aa 

+    [-(2+   ^  an  +   J  nAa)QQ   -  |  Qq   +   2   |]e 

-.         p  -nQ„Aa 

A2   =   ^  nQ^(4   -   e        "      ) 

1  P        7  -2nQp,Aa 

A^   =    [i  n(Aa-a)QQ   -   ;[^  Q^le 

1   ^2   -2nQoAa 
A5  =  •2  hQq  e 

It  can  be  shown  directly  from  (6.I6),  (6.I8)  and 
(6.20)  that 

A-,(x  =  Aa)  =  A2Cx=Aa)  , 


(6.22) 


A2(x=a-Aa)  =  Ao(x=a-Aa)  . 
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We  remark  that  for  the  hard  core  tail  potential  which 
we  postulate  over  a  <  |x|  <  a+Aa,  the  first  derivative  of 
A(x)j   0  <  X  <  a  J   does  not  have  to  be  continuous.  Indeed, 


(6.23) 


I  2   -nQn^^ 

A2(x=Aa  )  -  A-j_(x=Aa  )  =  -  2nQQ  e    ^ 


A-^(x=a-Aa  )  -  Ap(x=a-Aa  )  =^nQQ  e 


C. 


Calculation  of  Direct  Correlation  Function, 


We  are  interested  in  obtaining  b, (x)  over  the  "middle" 
interval,  Aa  <  x  <  a-Aa,  since  this  interval  is  dominantly 
larger  than  the  other  two  for  Aa  small,  and  also  the  comparison 
between  the  first-order  theory  and  exact  theory  will  be 
optimum. 

Referring  to  (5.53),  we  work  out  the  following 
Integrals,  for  Aa  <  x  <  a-Aa: 


(6.24a) 


A(C)  dC  =  vIq  +  y-j^x  +  v^x' 


(6.2ilb) 


(6.24c) 


X 
r 


?A(?)  dc  =  Pq  +  v^x^    +  v-^x^ 


0 
a-x 


A(?)  d?  =  e„  +  e-,x  +  c^x' 


(6.24d) 


a-x 

r 


—    —     —   2    —   3 
^A(?)  dc  =  e^  +  e-,x  +  e^x   +  e^x-^  , 


'0 
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in  which 

(6.25a) 

(6.25b) 

(6.25c) 

(6.25d) 

(6.25e) 

(6.25f) 

(6.25g) 


Pq  =  (Aa)-C^  +  I  (Aa)^-C2  "  ^^^^^l  "  I  (Aa)^B. 


y.  =  B- 


Vo  = 


2  ^2 


=  |(Aa)^C^  +  ^(Aa)^C2-  |(Aa)^B^-  ^i^a)^B^ 


yo  = 


U^  = 


e^  = 


|b. 


3  ^2 


(Aa)B^  +  |[2aAa-(Aa)^]B2 


+  [ 


-Ao      An      A^  1      ? 

-   — -  +   -^2_  (i+anQg)]-  (Aa) (A^+aA2)+  |(Aa)  A2 


'"^O    '^^O    n-Q 


0 


A,   -nQ„Aa    A,,    A^  nQ„Aa   ,.^s     nQ^^Aa 


'0 


'0    n  Q 


0 


'0 


(6.25h) 
(6.25i) 


£-,  =  -  B,  -  aB2 


r   =  1  R 
^2    2   2 


(6.25J)        Eq  =  |[2aAa-(Aa)^]B^+  j[3a^Aa  -3a(Aa)^+(Aa)^]B2 


+   -^   [A,(l-anQQ)+  A^(l+anQ^)+  ^(a^n^Q^   +  2anQn+  2)] 


2^2  ^"3 
n  Qq 


A 
'0'  ■  nQ 


0 


aAa(A^+aA2)+  |(Aa)^(A^+2aA2)-  ^(Aa)^A2   + 
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A,  -nQ.ia   ,.  ,    -nQ.Aa 

-,  A^   P  P  P  nQ„Aa 
P^  [A.d+anQ^)  +  -^(a^n^Q^  +  2anQ„+  2)]e   - 

(6.25k)        e^  =  -aB^  -  a^B2 

(6.251)        ^2  "  I  ^1  "^  ^^2 

(6  .  25m)        i";.  =  -  ^5-  Bp  . 

3    ^ 
We  now  construct  the  quantity  b-,  (x)  -  1   B„x   appearing 

^      0   ^ 
in  (5.53),  for   Aa  <  x  <  a-Aa.   Using  (6.24)  there  follows 

3        3 
(6.26)         b-,(x)  -  I   BgX^   =  I   V  x^  ,       Aa  <  X  <  a-Aa, 

in  which  the  B  are  given  by  (5-52)  and  (5.67),  and  the 
Vj  are  merely 

(6.27a)   Vq  =  -  B^-  2n^Q^yQ-  2nQQ(l+anQp)eQ  +  2n^Q^FQ 


(6.27b)   v-L  =  -  B2+2n^QQyQ+  2n^QQeQ-  2nQQ  (l+anQ^)  6-^+  2n^QQe^ 

2  2      ^2 

(6.27c)  v^   =    2n  Qq^^i"  ^12"^  ^1  "  ^^2  ^  ^2  ~  nQ~^ 

(6.27d)   ^^3  =  2  n^Q^iv^   -   "H^   +    e^   +    e^)    . 
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By  substituting  from  (6.25)  and  (6.19),  and 
manipulating, (6 . 27c )  and  (6.27d)  become,  respectively 

(6.28)  ^2  =  -  2n'^Q^(l  +  anQQ)(l  -  e    "   ) 

and 

i|   n  h       -nQ  Aa 

(6.29)  ^3=1  n  Qq^I  -  e      ) 


from  which 


(6.30) 


v^  3(l+anQQ) 


V3       2nQQ 


Inspection  of  (5.52)  enables  us  to  deduce  that 


V       B 
(6.31)  ^  -    ^ 


V3     B3 

But  (6.31)  is  the  necessary  condition  that  b-,  (x)  be  linear 
for  Aa  <  X  <  a-Aa   (refer  to  (6.26)),  and  strongly  indicates 
this  . 

Thus,  we  conclude  that  the  first-order  theory  here 
yields  a  straight  line, 

(6.32)  h-^(x)    =    (Vq+Bq)  +  (v-^+B-^)x  ,     Aa  <  X  <  a-Aa, 
in  accord  with  the  exact  theory,  viz., 

(6.33)  ^exact^^^  "  "  Q^^-^x)  ,        0  <  x  <  a+Aa 


with 


Q"-*"  =  1  -  n(a  +  Aa)  , 
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and  we  are  comparing  the  lines   b-,(x)  versus 

[C    .  (x)-  c„(x)]   over  Aa  <  x  <  a-Aa   (Aa  small  but 
exact      0 

finite)  . 

The  slopes  and  intercepts  of  b  (x)  and 
[c    (x)-c^(x)]   are  unequal,  for  finite  Aa,  but  the 
interesting  feature  is  that  the  first-order  perturbed 
equation  does  yield  the  correct  functional  dependence 
(linear)  over  the  interval  Aa  <  x  <  a-Aa  of  comparison 
with  exact  statistical  mechanics.   Even  though  this  might 
be  expected,  as  4>,  is  a  replica  of  a  potential  known  to 
provide  a  direct  correlation  function  which  is  linear, 
its   demonstration  forthwith  is  reassuring  and  lends 
credence  to  our  technique. 

In  the  limit  Aa  ^  0,  the  slopes  and  intercepts  of 

both  b-,(x)  and  [c    (x)  -  c„(x)]   approach  zero,  and  in 
-L  e  X  •        u 

t  , 

the  same  way.   That  is,  b.,  (Aa)  equals  c    (Aa)  as  the 

replica  sphere  vanishes. 
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VII.    Applications. 

A.     Hard  Cores  with  Surface  Adhesion. 

The  extended  hard-core  kernel,  A(x,x')j  of  equation 
(6.8)  can  be  used  to  obtain  the  first-order  change  in 
direct  correlation  function  for  the  sticky  hard-core  ("glue") 
model  of  Baxter  [l6].   In  this  model,  the  potential  consists 
of  a  hard  core  plus  a  rectangular  attractive  well  for  which 
a  certain  limit  is  taken  to  make  the  range  of  the  well 
become  zero  and  its  depth  infinite.   This  is  the  dual  of 
the  Kac  model  [17],  in  which  the  tail  potential  has  infinite 
range  and  vanishing  amplitude.   Both  models  exhibit  a 
first  order  phase  transition  when  solved  rigorously. 

In  the  glue  model, 

(7.1)  e[slue)(^,^  ^  A[6(x'-a)  +  6(x'+a)] 

where  A  is  a  normalization  constant.   From  (6.1)  and  (6.8), 


(7.2)  6c^^^^(x) 


A(x,x')  e|^^^(x')  dx'  . 


Thus 


(7.3)  6cS^"^(x)   =   2  A  A(x,a)  . 


Using  (6.7),  the  result  is 

(7.M    6cSl^e(x)  =  2AQn[-nc„(x)  +  ^  6(a-|x|)]  , 


with-  Qq  =  (1-an) 


-1 
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p"  "I  1  ]  f^ 

The  expression  (7.^)    for  6c^    (x)  Is  of  the  same 
type  as   (c^fj^f^(x)  -  c„(x)).   From  (2.21),  there  is 
the  equation  of  state: 


(7.5)   1^  3P  =  1  -  n  I  c(x)  dx  =  1^  BPq  -  n  J  6c^^'^^{x)    dx 


9  2 

where  -^  3Pf^   is  the  hard  core  contribution,  equal  to  Qq  . 

Substituting  5c^"'"^^(x)  from  (7.^)  into  (7.5),   we  obtain 


(7.6)  |_  gp  =  (1  _  2nAQo)  1^  BPo  , 


with  a  van  der  Waals  loop. 

The  use  of  A(x,x')  in  (7.2)  to  obtain  6c^'^^®(x)  is 
not  merely  ad  hoc;  this  technique  can  also  be  applied  to 
other  short-range  tall  potentials.   It  is  the  restriction, 
an  £  1 ,  on  the  size  a  and  density  n  of  the  cores  which 
determines  the  scope  of  the  method;  the  strong  singularity 
of  A(x,x')  at  |x'|=  1/n   corresponds  to  close  packing  and 
consequently  imposes  a  limitation  on  the  range  of  the  tail 
potentials  which  can  be  treated. 
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B.      The  Debye-Huckel  Model  with  Hard  Core  Insertion. 

In  the  previous  chapter  we  established  that  b, (x) 
is  linear  for  an  extended  hard  core  of  finite  range.   Since 
the  theory  is  linear  in  e, ,  scaling  arguments  indicate  that 
soft-core   tail  potentials  also  show  a  linear  region  for 
b-,  (x)  ,  over  0  <  x  <  a.   Since  arbitrary  short-range  (})-,(x)'s 
can  be  composed  approximately  of  incremental  cores,  it  is 
expected  that  b-,  (x)  will  remain  linear  for  such  (j)-,(x)'s. 
For  long  range  tails,  however,  such  arguments  cannot  be 
sustained  (the  condition  an  £  1  would  be  violated). 

However,  if  the  phenomenon  of  linearity  is  extrapolated 
as  a  model  for  arbitrary  tail  potential,  and  the  P.Y. 
equation  is  further  simplified  to  the  Debye-Huckel  model  with 
hard  core  insertion  (extended  mean  spherical  model  [l8,19]), 
a  number  of  explicit  solutions  can  be  obtained. 

The  D.H.  model  with  hard  core  Insertion  is  defined  by 
the  following  equations: 

(7.7a)         g(x)  -  1  -  c(x)  =  n(g-l)  *  c  ,   -oo  <  x  <  «> 

(7.7b)         c(x)  =  -3  ())(x)  ,  |x|  >  a 

(7.7c)         g(x)  =0       ,  |x|  <  a   . 

The  existence  of  a  hard  core  is  indicated  by  (7.7c). 
This  hybrid  model  is  computationally  simpler  than  the  full 
P.Y.  equation,  and  is  useful  now  since  one  of  the  two  unknowns. 
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namely,  c(x)  for  |x|  <  a  (the  other  unknown  Is  g(x),  for 
|x|  >  a)   can  be  assumed  linear  as  a  result  of  the  prior 
approximation  that  c(x)  in  vicinity  of  origin  is  shielded 
from  the  explicit  tail,  and  feels  only  the  net  pressure. 
Thus,  we  write 


(7.8) 


c(x)  =  a  +  Y  X  > 


0  <  X  <  a  , 


in  which  the  constants  a  and  y  are  to  be  determined  by  the 
continuity  conditions  on  c  and  g  in  this  model.   Thus, 
from  (7.7a,c)  for  y  >  0: 

oo  oo 

f 


(7.9)   g(y)  -  c(y)  =  (1  -  2n  J  c(x)  dx)  +  n 

0 


+  n 


c(x-y)g(x)dx  +  n 


c(x+y)  g(x)  dx 
c(x-y)g(x)dx, 


a+y 


a+y 


Differentiating  (7.9),  there  follows 

00  00 

(7.10)  g'(y)-c'(y)  =  n  J  c ' (x+y )g(x)dx  -  n  J  c ' (x-y )g(x)dx 


-nc(a  )g(a+y)-  n 


a 


a+y 


a+y 


c ' (x-y)g(x)dx  +  nc(a  )g(a+y)  , 


and 


(7.11)  g"(y)-c"(y)  =  n    c"(x+y)g(x)dx  +  n   |  c" (x-y )g(x )dx 


a+y 


+nc'(a  )g(a+y)-  nc(a  )g'(a+y)+  n 


a+y 


c"(x-y)g(x)  dx 


a 


-  nc'(a  )g(a+y)  +  n  c(a  )  g'(a+y) 
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Thus  J  as  y  approaches  zero,  (7.9)-(7.11)  take  the  form 


00 

r 


(7.12) 


-c(0  )  =  (1  -  2n 


c(x)  dx)  +  2n 


oo 

r 


c(x)g(x)  dx 


(7.13)         -c'(0"^)=  n(c(a")  -c(a"^)  )g(a"^) 


(7.14)         -c"(0^)=  n(c(a")  -  c  (a"^)  )g' (a"^) 


+        +  r 

-  n(c'(a  )  -  c'(a  ))g(a  )  +  2n   c"(x)  g(x)  dx  , 


respectively.   Since   g-c   and  Its  first  derivative  are 
continuous  at  a,  and  g(a~)  =  g'(a~)  =  0,  we  have 


(7.15) 


g(a  )  =  c(a  )  -  c(a  ) 
g'(a^)  =  c'(a"^)-  c'(a") 


Using  (7.15),  (7.12)-(7.l4)  simplify  to 


(7.16) 


(7.17) 


,+ 


-c(0  )  =  (1  -  2n 


c(x)  dx)  +  2n   c(x)  g(x)  dx 


_c'(o'^)=   -  n  (c(a'^)  -  c(a"))' 


(7.18) 


■c"(0'^)  = 


2n 


c"(x)  g(x)  dx 


It  is  noted  that  for  no  tail  at  all,  <))  =  0 ,  (7 .  l6)-(7 .  l8) 
correctly  imply 
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(7.19a) 


(7.19b) 


•Cq(0'^)   =   1  -  2n  I  Cq(x)  dx  =  Qq 


0 


■Cq(0'^)   =   -  "  Qq 


(7.19c) 


-0^(0"^)   =  0   . 


Long-Range  Forces. 


We  now  carry  out  an  explicit  solution  for  an  attractive 
Yukawa  tail  potential. 
The  potential  is 


(7.20) 


<t)(x) 


=   -  A  e-^l^l 


for  |x|  >  a. 


From  (7.7b), 


(7.20') 


c(x)  =  BA  e 


-Ax 


X   >  a, 


Clearly, 


(7.21) 


c"(x)  -  A^c(x)  =  0   , 


X  >  a, 


It  follows  from  (7.8)  that 


(7.22)     c(0)  =  a  ,   c'(0)  =  Y  ,   c"(0)  =  0 


By  means  of  (7.l8),  (7.21)  and  (7.22)  we  find 


(7.23) 


c"(x)  g(x)  dx  =  0  ;       c(x)  g(x)  dx  =  0 


Using  (7.8)  and  (7.20')  indicates  that 
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(7.24) 


c(a~)  =  a  +  ay  ,    c(a  )  =  3A  e     , 


and 


(7.25) 


,    s    ^  ,1   2  ^  BA   -Xa 

c(x)  dx  =  aa  +  p-  ya  +  — r—  e 


0 


Substituting  (7.22),  (7.24)  and  (7.25)  into  (7.16)  and  (7.17) 
yields  two  equations  for  a  and  y  • 


(7.26a) 


-  a  =  1  -  2n(aa  +  p"  Y^   +  — r—  e    ) 


(7.26b) 


y  =  n(a  +  ay  -  3A  e~   ) 


Referring  to  (2.21),  (2.22),  (7.l6),  (7.22)  and  (7.23), 
there  exists  the  relation 


(7.27)    1^  (6P)  = 


n  ^T 


=  1  -  2n 


c(x)  dx  =  -  c(0)  =  -  a 


Equations  (7.26a,b)  can  be  solved  readily.   First,  from 
(7 . 26a)  ,  we  have 


(7.28) 


ay  = 


1    2  3A  -\a.    _L   /  1    o^ 
—  -  — r—  e    )  +  a(— -—  -  2) 


an    aX 


an 


Substituting  (7.28)  into  (7.26b)  in  order  to  eliminate 
y  leads  to  the  following  quadratic  equation  for  a  : 


(7.29) 
with 


pa   +  qa  +  r  =  0 
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p  =  (1  -  an) 
(7.30)         q  =  1  +  ^^  (an  -  l)(aX  +  2)6"^^ 


The  two  roots  of  (7-29)  are 

1  +  2n3A  (an-l)(aX+2)e-^^ 

(7.31)   a  =  - p 

2(l-an)'^ 


1/2 
(1  _  ilJlM  (i-an)[aX(2an-l)+2an]e"^^} 

2(l-an)^ 


Since  a  must  be  real,  we  require 


(7.32)         ^i^  (l-an)[aX(2an-l)  +  2an]  e"-^^  <  1 


The  sign  in  front  of  the  square  root  in  (7-31)  must  be 
chosen  so  that  the  compressibility  is  positive,  km  >  0. 

By  placing  A  =  A  in  the  attractive  Yukawa  potential, 
equation  (7.20),  and  then  taking  the  limit  X  ->  0 ,  we  recover 
the  tail  potential  of  the  Kac  model  [17],  which  is  known 
to  lead  to  a  phase  change  described  by  the  van  der  Waals 
equation: 

(7.33)         6P  =  Y^  -   Bn^  , 
from  which 
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(7.34)         f-  (3P)  =  ^  -  ^n3(l-an)^ 
^"^  (1-an)^ 


Evaluating  (7.31)  in  the  limit  X  -^  0 ,  we  find 

1  +  4nB(an-l)  +  /l-8an^3(l-an) 


(7.35)  I-   (3P)  =  -  a  = 


^"^  2(l-an)2 


At  very  high  temperature,  the  properties  are  determined 
by  the  cores.   In  the  limit,  3-^0,  (7.34)  and  (7.35)  agree. 


(7.36)         lim  |-  (3P)   =  qI    . 


For  finite  3,   >  3q,   the  van  der  Waals  description  is 

recovered . 

The  above  technique  obviously  applies  to  double 

Yukawa  tails.   Such  potentials,  ()) ,  which  in  three 

dimensions  have  the  form 

-,      ~^-\^  -X^r 

(7.37)         (J)(r)  =   ^  (A^e   ^  -   k^   e      ^    )    ,      0   <   v   <   -> , 

simulate  intermolecular  potentials  quite  well. 
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Appendix  A.   Evaluation  of  rCx)  from  r"(x). 

1.    Expression  of  fCx)  in  Terms  of  a  Single  Constant. 

Prom  (5.36)  and  (5.39),  for  -a  <  x  <  2a, 

(A.l)   r"(x)  =  6"(x)  -  2n^QQ6(x)  -  2n^QQ6(x-a)  +  2nQQ6'(x-a), 

Writing  c„(x)  (refer  to  (5.11))  as  a  continuous 
function  plus  a  step  function  leads  to  the  conclusion  that 
Cq(x)  can  be  expressed  in  terms  of  delta-functions, 

(A. 2)   Cq(x)  =  2nQQ6(x)  -  nQQ6(x-a)  +  Qq6 ' (x-a)  . 

(A. 2)  enables  (A.l)  to  be  written 

(A. 3)   r"(x)  =  6"(x)  -  6n^QQ6(x)  +  2ncQ(x) 

or 

(A. 4)   (r(x)  -  6(x)  -  2ncQ(x))"  =  -  6n^Q^&{x)    . 

Since  the  left-hand  side  of  (A. 4)  is  an  even  function  of  x, 
there  must  be  no  odd  functions  appearing  in  the  second  anti- 
derivative.   Thus,  a  single  (finite)  constant,  denoted  by  K, 
appears  in  the  result 

(A. 5)   r(x)  =  6(x)  +  2ncQ(x)  -  3n^QQ | x |  +  K  ,   (-a  <  x  <  2a), 
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Calculation  of  the  Constant 


From  (3.53)  and  (5.2),  there  follows 


(A. 6)   (r(x)  -  6(x))  «  (6(x)  -  ncQ(x))  =  n(gQ(x)-l)  +  ncQ(x), 

valid  for  all  x.   In  particular,  at  x  =  0 

a 
(A. 7)   (r(0)  -  6(0))  -  n 


CQ(x')(r(x')  -  6(x'))  dx- 


-a 


=  -  n  +  nc„ (0)  . 

Using  (A. 5)  to  introduce  K  into  (A. 7)  leads  to  the  relation 

a 

-3n2 


(A. 8)   K  =  g-  k^[-  n  -  ncQ(O)  -  Sn^Qg 


x  c^(x)  dx 


0 


+  I4n' 


? 


0 


Cq(x)  dx] 


in  which  the  isothermal  compressibility,  k^,  occurs  due  to 
(2.21).  In  view  of  (5.11),  all  terms  in  (A. 8)  are  readily 
calculated.   The  result  is 


(A. 9)   K  =  a  n^Q^  (6  -  6an  +  i^  a^n^  -  a^n^)  , 


n  , 
since  |-  krp 


,-2 


for  hard  cores 
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Appendix  B.   Calculation  of  r_2(x-2a). 


From  (5.33),  (5.36),  and  (5.37),  the  only  non-zero 

ti       II 
terms  in  r"(x)  for  2a  <  x  <  3a  are  r^Cx),  r_-|_(x-a),  and 

r_2(x-2a) . 

Referring  to  equation  (A. 3),  we  have 

(B.l)   r"(x)  =  6"(x)  -  6n^QQ  6(x)  +  2ncQ(x)  +  r"2(x-2a), 

(2a  <  X  <  3a) 


The  term  r  p(x-2a)  is  the  additional  piece,  addended 
when  X  >  2a.   From  (B.l),  there  follows 


(B.2)   r(x)  =  6(x)  +  2nCQ(x)  -  3n^Q^ | x |  +  K 


r_2(x-2a)  d^x  , 


(2a  <  X  <  3a) , 


with 


(B.3)   r_2(x-2a)   =  II  r"2(x-2a)  d^x 


From  (5.37), 


00             /   2  2    2x 
(B.il)   r  2(x-2a)  =  -  ^i    g-ik(x-2a)  0 


dk 


Using  a  partial  fraction  decomposition,  we  find 


-7^- 


1,           ■:,':,                       6n^Q^ 
(B.5)   r_2(x-2a)  =  3n^QQ6(x-2a)  ^ 


-lk(x-2a)    dk 
^  (nQQ-lk) 


4  4  00 
+  ^   ^Q   f  g-ik(x-2a)  dk_ 


27r 


(nQg-ik) 


,  (x  >  2a) . 


The  Integrals  are  readily  evaluated  by  contour  techniques. 

(B.5)  can  be  written 

„  P  P  ^    ^      -nQ„(x-2a) 

(B.6)   r_2(x-2a)  =  3n'^QQ6(x-2a)  -  6n^Q^  e   ^ 


4  4        -nQQ(x-2a) 
+  n^Qo(x-2a)  e   ^ 


(x  >  2a) 


Prom  (B.3)  and  (B.6),  we  obtain  the  result 


=  3  ^2^2 


(B.7)   r_2(x-2a)  =  ^  n''QQ|x-2a|  -  4n  Qq  e 


-nQQ(x-2a) 


P  P        -nQ„(x-2a) 
+  n'^QQ(x-2a)  e   ^ 


(x  >  2a) . 
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